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RESULTS FROM THE NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION REMOTE SENSING EXPERIMENTS 
IN THE! NEW YORK BIGHT - APRIL 1975 
John B. Hall 5 Jr. and Albin 0. Pearson 

ABSTRACT 

A cooperative operation he ween the National Aeronautics and Space 
Administration (NASA) and the National Oceanic and Atmospheric Administration 
(NOAA) ■was conducted in the New York Bight during April 7-17 } 1975* to eval- 
uate the role of NASA remote sensing technology to monitor ocean dumping. Six 
NASA remote sensing experiments were flown on the C-5U, U~2 s and C-130 NASA, 
aircraft, while NOAA obtained concurrent sea truth information using heli- 
copters and surface platforms. The experiments included l) a. Radiometer/ 
Scatterometer (rADSCAT), 2) an Ocean Color Scanner (OCS), 3) a Multichannel 
Ocean Color Sensor (MOCS), four Easselblad cameras, 9) Ebert 
spectrometer, and 6) a Reconafax IV infrared scanner and a Precision 
Radiation Thermometer (PRT-5)* 

The purpose of this report is to present the results of these experiments 
relative to the use of remote sensors to detect, quantify, and determine 
the dispersion of pollutants dumped into the New York Bight. Descriptions of 
the experiments, experimental methods, data analyses techniques , end 
significant results for each experiment are given in separate sections of 
the report . 
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RESULTS FROM THE NATIONAL 
AERONAUTICS AND SPACE ADMINISTRATE ON 
REMOTE SENSING EXPERIMENTS IN 
THE NEW YORK. BIGHT - APRIL 7-17, 1975 

COMPLIED BY 

John B« Hall Jr* and Albin 0. Pearson 
SUMMARY 

Six remote sensing experiments were conducted "by the National Aeronautics 
and Space Administration (NASA) in con junction with the National Oceanic and 
Atmospheric Administration (NOAA) in the New York Bight between April 7-17 s 
19T5j to evaluate the role of NASA remote sensing technology to monitoring 
ocean dumping. Sixteen remote sensors were flown on the C-5**, U-2, and C-130 
NASA aircraft , while NOAA obtained concurrent in situ sea truth data using- 
helicopters and surface platforms. The six primary sensor experiments 
included (l) a Radiometer Seatterometer (RADSCAT) , (2) an Ocean Color Scanner 
(OCS), (3) a Multichannel. Ocean Color Sensor (MOCS), (4) four Hasselhlad 
cameras, ( 5 ) an Ebert spectrometer, (6) a Reconofax IV infrared scanner and a 
Precision Radiation Thermometer (PRT-5)* 

The purpose of this report is to present the results of these 
experiments. Descriptions of the primary sensors, experiments! methods, 
and data analysis . techniques along with significant results for each 
experiment, are given in separate sections of this report. 


INTRODUCTION 


The National Oceanic and Atmospheric Administration (NOAA) 4 * in 197^ 
initiated a 7-year Marine Ecosystem Analysis (MESA) Program to obtain environ- 
mental information in selected marine locations* The New York Bight is the 
first site selected for intensive study under the MESA program. The program 
is focused to provide information to better understand* utilize, and minimize 
man’s impact on the New York Bight. A more efficient utilisation of the 
Bight resources requires techniques for a rapid and accurate assessment of 
the effects of man’s activities on the Bight ecology. This requirement, 
coupled with the dynamic nature of the marine environment, accents the need 
for the synoptic spatial and temporal advantages of remote sensing systems. 
Accordingly, as a cooperative effort in the MESA program, NQAA requested that 
NASA evaluate the use of remote sensing technology to define circulation 
features in the New York Bight and to apply this technology for monitoring 
and managing ocean dumping* The first step of this cooperative effort was 
accomplished in the New York Bight during April 7-17, 1975* The area 
selected by NOAA for this joint NOAA/NASA evaluation was primarily the apex 
of the New York Bight as shown in figure 1. The New York Bight extends from 
Cape May, New Jersey, to Montauk Point, New York, and seaward to the. edge 
of the continental shelf (200 meter depth). The apex is bounded on the north 
by Long Island, on the south by latitude JiO^O’N, and on the east by longitude 
meridian 73°30’W. Presently, within the apex, sewage sludge is dumped at a 

*The authors are grateful to John W. Sherman and James B. Zaitzeff of NOAA 
for their guidance and coordination of NOAA activities in implementing this 
effort . 
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location about 18 km (9,7 nautical miles (n. mi.)} south, of Long Island. Acid 
■wastes are dumped at a location about 10 km ( 5 *^ n. mi.) southeast of the sewage 
sludge dump . site. Six remote sensing experiments were conducted by MSA 
using l 6 remote sensors flown onboard aircraft platforms to obtain sea surface 
information relative to these dump sites. Concurrently NOAA obtained in situ 
"sea truth 11 information with a combination of helicopters and surface plat- 
forms. NASA operations associated with these experiments are given in 
reference 1 . 

The purpose of this report is to present the results of the NASA remote 
sensing experiments. The following sections of this report contain 
descriptions of the remote sensors , experimental methods , data analysis 
techniques, and pertinent experiment results. 


Reference 


Usry, J. W. and Hall, J. B. 5 Jr.: National Aeronautics and Space 

Administration Operations Remote Sensing Experiments in the Nev York 
Eight, April 7-17* 1975. NASA TMX-72802, November 1975- 






Figure 1. - MESA Defined New York Bight and Apex. 
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QCEM COLOE SOAUIfER (OCS) 
BY 


WABEEN A. HOVIS 


OCEAN COLOR SCANNER (OCS) 
by 


Warren A. Hot is 
SUMMARY 

The MSA U-2 aircraft participated in the MESA experiments in the Hew 
York Bight during April 7-17, 1975- The OCS and four Mitchell-Vinten 70 -mm 
cameras were flora onboard the MSA U-2 aircraft on April 9, 13, and it, 
1975* The purpose of this report is to present the results of this 
investigation. 

EXPERIMENT DESCRIPTION 

The Ocean Color Scanner (fig. 1) ms flora on the U-2 aircraft at an 
altitude of 19.8 km (65,000 ft) and ms supported by four 70-mm Mitchell- 
Vinten cameras. The OCS is a 10-channel multispectral scanner ranging from 
^33 to 772 nanometer (nm) center wavelength with a 90° total scan angle 
and spatial resolution of 3-5 milliradians. Spectral bands and saturation 
■radiance values used in these investigations are shown in table 1. The 
radiance for saturation shown in table 1 is for a gain of one for each 
channel. Gain may be increased separately for each channel in steps of 1.5, 
2, and 3 to allow for changing Sun angle due to seasonal changes or time of 
flight. All channels, except 9, are optimized for water scenes, including 
atmospheric .backs catter as seen from an altitude of 19.8 km. Channel 9 has 
the same level of gain as the Landsat band 6 that covers 700 to 800 nm. 
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Four 70-mm Mitchell-Vinten cameras were flown with the OCS to provide 
high resolution photography to aid in locating the U-2 position. Three of 
these cameras contained hlack and white film, and one contained color film. 

A photograph of this system is shown in figure 2. Specific cemera details 
pertinent to this experiment are shown in table 2. 

DATA ACQUISITION 

Flight lines for the U-2 are given in reference 1. These lines were 
established to allow the aircraft to fly directly toward, or away from the 
Sun, to avoid Sun glint in the scanner imagery. The morning flights were 
in a general southeast to northwest direction while afternoon, flights were 
generally southwest to northeast. A slight "fan" effect is seen due to 
allowance for Sun motion during the time of each flight sequence. The actual 
flight lines flown, flight altitudes , flight line locations, and sensor 
operation times are given in reference 1. 

DATA ANALYSIS AND RESULTS 

The flight on April 9 was successful for U-2 operations. Surface truth, 
however, was not available due to lack of communication between the ship 
and helicopter operations in the Bight, and the U-2 base at NASA Wallops 
Flight Center. The U-2 decision to fly was based upon the weather forecast 
available at Wallops. 

The photographic and scanner data were both of excellent quality for 
this flight. The scanner data were processed to a pictorial product using 
two different enhancement techniques and, despite the lack of surface truth, 
still show some interesting information as shown in figure 3. Figure 3 is a 
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n conventionally 1 ' enhanced image. It -was made from art enhancement where a 
direct current offset is applied to each of the channels used in making the 
image to allow for atmospheric backscatter. The amount of offset is 
wavelength dependent. The image in figure 3 w .s made from the enhanced data 
from channels 465, 582 s and 662 nm. The image shows the Slew Jersey shore 
from Barnegat Inlet to just south of Sandy Hook Point and a portion of the 
south shore of Long Island. A new acid dump Is seen as a sharp n U n , while 
the older dumps are shown drifting to the west. The westward drift is 
characteristic of all imagery obtained of this area. In addition, the older 
dumps are shown as being diffused in a southerly direction. Also shown In 
figure 3 is a square-shaped dump to the north of the acid waste dumps. This 
feature is a sewage sludge dump that was observed in its early stages 
during the morning overpass at approximately 1450 g.m.t. The brightest part 
of this feature, the north edge in the northwest corner, is the Initial 
dumping area. This indicates that the dump is not uniform with time and that 
the concentration is higher in the initial dump area or that the sludge is 
not uniformly mixed. It appears that the material dumped first Is brighter 
in the spectral regions sensed than that dumped later. The bright area along 
the Hew Jersey shore Is the Hudson River plume. This is the normal flow 
pattern of the plume as seen in a number of images of this area. The high 
relative reflectivity of the plume indicates a high particulate loading which 
was confirmed by sea truth measurements on April 13. The plume appears 
rather featureless in figure 3. Further enhancement was performed by squaring 
all of the radiance values in the first enhanced tape. The image that 
results from this enhancement is shown in figure 4. In this image, 
gradations in the plume are more easily observed. An apparent layering on the 


east edge in the plume is observed as the plume bends after leaving the 
Hudson River. The April 13 flight ■was performed in the morning. Cloud cover 
precluded an afternoon flight. The position of the helicopter sampling grid 
was located on the scanner image, and radiance levels at most of the 
helicopter stations were determined. The radiance, in all channels, showed 
a general increase with increased sediment concentration. * Figure 5 shows the 
measured upwelling radiance as measured by the sc ann er versus the log of the 
measured sediment concentration. There is no immediate explanation for the 
four points at the high sediment levels that appear offset from the others. 

In an attempt to extract qualitative information from the radiometric data, 
the method of characteristic vector signatures was tried with the addition 
of variation from mean radiance as an attempt to compensate for atmospheric 
effects. Briefly, a mean overall radiance level is determined for the water 
in the scene, and the departure from the mean determined for each spectral 
band. As show in the formula in figure 6, the total radiance consists of 
the mean radiance, the atmospheric radiance, and a contribution due to 
sediment and chlorophyll. The vector signature can be either positive or 
negative which denotes either scattering brighter than the mean or absorption 
darker than the mean. The sediment signature shown at the top of figure 6 
is positive at all wavelengths which indicates no specific absorber, i.e, , a 
near grey color. The chlorophyll signature shows that some chlorophyll 
absorption can be detected in the plume even though chlorc£jhyll is a minor 
constituent. There is a reasonable degree of correlation between the total 

*The author is grateful to Dennis S. Clark and Terry A. If el son of IfOAA for 
the use of their unpublished sea truth information. 
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sediment and the chlorophyll concentration (fig. 7)« This is to he 
expected since the plume is ‘virtually all of organic origin. The chlorophyll 
content is measured in pg/& while the sediment content is measured in, mg/Ji. 

An attempt was made to relate measured radiance to chlorophyll content in the 
plume with only modest success as shown in figure 8. This result is expected 
when sea truth information shows that only a small portion of the plume is 
composed of chlorophyll. Chlorophyll has its specific absorption, while the 
major components of the plume hare spectral character. Correlation of the 
calculated and measured sediment concentrations is given in figure 9- Here 
the agreement is good and would give considerable confidence that a contour 
map of total sediment could be made by applying the sediment vector 
signature to the entire area of the plume. It is apparent from the plume 
analysis that it will be difficult to determine concentrations of minor 
constituents of a plume from radiometric data only. It may be possible, 
however, to infer the minor constituent concentration if its quantitative 
relation to the total sediment load is known. 

The magnetic tape recorder failed during the April lU flight. Failure 
was due to damage sustained in a previous shipment from Wallops Flight Center 
to the Ames Research Center. The photographic quality was good, but could 
not be used for quantitative analysis. 


Reference 

1. Usry, J. W. and Hall, J. B. , Jr,: National Aeronautics and Space 

Administration Operations - Remote Sensing Experiments in the New York 
Bight, April 7-17, 1975* MSA TMX-72802, November 1975. 









TABLE X.- OCEAN COLOR SCANTIER 


Channel 

Cent er 
Wavelength, 
(nm) 

Bandwidth 

(nm) 

Radiance for 
Saturation 
(Gain X 1) 

mv/ cm^ 

1 

433 

22.5 

4o.io 

2 

471 

21.5 

25.00 

3 

509 

27.0 

23.60 

4 

547 

24.5 

14.70 

5 

583 

25-0 

11.80 

6 

620 

26.0 

10.00 

7 

662 

22.0 

7.55 

8 

698 

20.5 

5.00 

9 

T33 

22.5 

11.90 

... 10 

,772 

23.0 

3-47- 




TABLE 2.~ MITCHELL-VIUTEIf CAMERAS 


Camera 

Number 

Lens 
Focal 
Length, 
cm(in. ) 

Film 

Format*. Number 

cm(in) 

Shutter 

Speed, 

sec 

Filter 

Number 

Spectral 

Hange 

Nanometer 

F 

Stop 

Spatial 
Resolution 
m (ft) 

1 ■ 

4.43 

(1.75) 

5.72 x 5.56 
(2 7 X 2 16 > 

Pari atomic 
X, 3400 

1/250 

Wratten 

12 

510-700 

9.6 

9.1 to 15.2 
(30 to 50) 

2 

4.43 

(1.75) 

5.72 X 5-56 

(e i z2 V 

Panatomic 
X, 3400 

1/250 

Schott GG 

475 + 

Schott BG 
18 

475-575 

5.6 

9.1 to 15.2 
(30 to 50 ) 

3 : 

4.43 

(1.75) 

5.72 x 5.56 

< 2 1 x 2 i6> 

Panatomic 
X, 3400 

1/250 

Scliott OG 
570 +• 
Schott BG 
38 

580 - 680 . 

5.6 

9.1 to 15.2 
(30 to 50) 

4 

4.43 

(1.75) 

5.72 X 5.58 

( 2 i X 2 k } 

Aerial 

Color 

50-242 

1/250 

None 

400-700 

3.5 

9.1 to 15.2 
(30 to 50) 
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SEWAGE SLUDGE 







CHARACTERISTIC VECTORS 



S=MEASURED RADIANCE <S> “SEDIMENT CONCENTRATION 

So = MEAN RADIANCE <C> = CHLOROPHYLL CONCENTRATION 
Figure 6 -CHARACTERISTIC VECTORS FOR SEDIMENT AND CHLOROPHYLL. 


MEASURED SEDIMENT (S) 



°° / 

*5 

1 

LOG CHLOROPHYLL CC) = 1.384 + 0,362 x LOG SEDIMENT CS) 

/ 

:i 

t 

i 


/ o 

< 

o 

/ 

1 

f 

o 

/ 

1 

f 

o / 

0 0 

J 

i 

i 

i 

i 

0 / 


o / 


t 

f 

0 o / 

HUDSON RIVER PLUME 

j 

> 

* 


APRIL 13, 1975 



LOG MEASURED CHLOROPHYLL CC) 

Figure 7 - LOG MEASURED CHLOROPHYLL VS LOG MEASURED SEDIMENT 






Preliminary Analysis of Ocean Color Spectra 
Measured Fit h an Ebert Spectrometer During the 
FASA/NOAA Spring 1975 Experiment in the Few York 3ight 

By 

James L. Mueller 
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PRELIMINARY ANALYSIS OF OCEAN COLOR SPECTRA 
MEASURED TOTH AN EBERT SPECTROMETER DURING THE NAS A/NO AA 
APRIL 1975 MESA EXPERIMENT IN THE NEW YORK BIGHT 

by 

James L. Mueller 


INTRODUCTION 

During the April 19T5 NASA/NOAA MESA experiment in the New York Eight, 

ocean color spectra were measured from a low flying aircraft. The short 

duration of availahle sampling flights precluded statistically correlating 

these ocean color spectra with simultaneous ground truth oh ser vat ions. 

Nevertheless, the ocean color spectra are suitable for comparison, using 

1 2 

characteristic vector analysis, with earlier results s and for studying the 
spectral signature of the acid waste pollutant routinely dumped in the New 
York Bight Apex. 

In radiometry, a ’’spectrum 11 is the distribution of radiant energy, or of 
some optical property of the medium, as a function of wavelength. The ocean 
color spectrum is defined as 


N. 

U = ,3 « 1, 2, . . ., 55, Cl) 

J N 


where is the upwelled radiance from the sea in the j 5 nanometer (mu) 
wide wavelength hand between 420 and 695 am, and the normalization factor N 
is the average radiance, i,e« 


55 

N = 1/55 2 N . 
Z=1 J 
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The spectrum U. must tie numerically parameterized "before it can be 

U 

efficiently correlated with other measurable properties of sea water. . 

] 2 

Mueller * 9 reviewed previously proposed parameter izations of ocean color 

3-6 

spectra and applied the method of characteristic vectors to the problem. 

It was shown for a sample of ocean color spectra measured off Oregon, that 

the first four characteristic vectors of the sample covariance matrix 

yielded a highly accurate parameterization. To compare Mueller’s earlier 

work with the MISA ocean color data, certain of those results must be recast 

in terms of the spectrum U. rather than in terms of reflectance 

J 


N. E 



where and E are the incident irradiance spectrum and its wavelength 

average. This change must be made because there is no consistent means of 
estimating E^ for the independent ocean color samples that will be 
examined. 

In the present report the characteristic vectors of U. spectra from 

t3 

the MESA experiment will be compared with similar results from two other 
independent samples of ocean color spectra measured elsewhere. The purpose 
of this comparison is to test the working hypothesis that ocean color spectra 
may be parameterized, globally, with a few standard, characteristic basis 
vectors. 


In addition, the spectral signature of the acid waste material routinely 
dumped in the Hew York Bight will be examined* It will be demonstrated that 
this material has a spectral signature that is uniquely different from any 
other ocean color spectrum yet observed. 




CHARACTERISTIC VECTOR ANALYSIS OF COLOR SPECTRA 


The use of characteristic vectors to parametrise ocean color spectra 'was 

12 8 

recently discussed by Mueller * . Morrison gives a comprehensive discussion 

of the method within the general context of multivariate statistical analysis. 

th 

A. spectrum U , ik® n of a sample of II ocean color spectra, is 
parametrised by its first K principal components 


Y, = 2 e T .{U. - M.); k = 1, 2, . 

kn kj v jn j * 


. K & n = 1, 2j . , . 3 H s 


( 2 ) 


where e is the characteristic vector of the sample covariance matrix, 
M. is the sample mean spectrum, and the index j denotes one of 55 discrete 

t3 

wavelength bands. The inverse approximation relating the spectrum to 


its first K principal components Y^ n 5 s 


K 

U - M + £ y e ; 3 - 1, 2, 

J n <J kj 


, J & n = 1 , 2 , . . . , N . 


(3) 


The characteristic vectors e, . are extracted from the sample' covariance 

kj 

matrix such that e_ . is alined in the direction of maximum sample 

tJ 

dispersion (i.e., it accounts for the largest possible proportion of sample 

variance). Then, e^. is extracted in alinement with the largest remaining 

3 

proportion of dispersion (sample variance) orthogonal to e p j * -Aiui ea °k 
succeeding unique characteristic vector accounts for a smaller percentage of 
sample variance than any of its predecessors. Tile accuracy with which a 


spectrum U. (j = 1, 2, . . 55) is parametrized by its first K principal 

jn 

components (k - 1, . . . , K), therefore, is measured by the cumulative 
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percentage of sample variance accounted for "by the first K characteristic 

vectors e. . . 

kj 

The t, Y. percent significant degrees of freedom 51 in ocean color spectra of 
a given sample will-. "be defined as the independent scalar parameters which 
retain a specified cumulative percent (y) of total sample variance* For each 
of the three samples considered here, there are five 99 percent significant 
degrees of freedom contained in ET (overall brightness) and the first four 
principal components (Y^, Yg, Y^, Y^) (spectral shape parameters). 

Thus far characteristic vectors applied only within the context of a 
single sample of ocean color spectra have been discussed. New some elementary 
relationships between characteristic vectors of two independent samples of 
ocean color spectra must he considered. For convenience, consideration will he 
given only to samples with five 99 percent significant desgrees of freedom. 

The characteristic vectors determined front a sample ot of ocean color 
spectra are related to those determined from sample 3 by the equation: 



t c“ 
k=i klB 


(3 , o. 

kj xj|S 


ih) 




.a 


where e. . D is an error vector describing an orthogonal discrepancy 
ijp 


between the two vector systems, and 


55 


C ki(3 e ij * e kj k " 1] 2) 3) k) • 

t) -L . 


(5) 


The mean ocean color spectra of the two samples are related by the equations 


iC “ MT + 2 b“„ eL + p a 


■j d * Ttp ^ 


( 6 ) 
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where p® a is an error vector similar to z\ 


a 


dB 




and 


,a 


B ^ = k = 1 > a * 3 * 4 ' 


(T) 


It follows directly from equations (2) and (4) through (?) that the principal 
components of an ocean color spectrum are related, in the two characteristic 
vector systems, hy the equation 


Y? = l c{\ R + 2 C® B<* -fr <5? 
i ± kip k k=:L kiP kp iP 


( 8 ) 


■where 6. Q are residual errors due to any unresolved orthogonality "between 
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the characteristic vectors of the two samples. 


When the spectra of sample a are parametrised with the sample P 


characteristic vectors, the cumulative proportion of total sample variance 
retained "by the cross-sample parametrization is given "by 


4 4 2 

V° = Z r“ [ I (c“ e ) ] 

3 1=1 1 k=l 1513 


(9) 


where v? is the proportion of sample variance associated with the. i 


principal component I, , i.e. , as determined with its own a-sample 
characteristic vectors. 


The method of characteristic vectors may also he extended to study the 


spectral signatures of specific pollutants as contained in a spectrum of U^. 


This is possible, however, only if color spectra of water containing that 
pollutant were not included in the sample used to derive the characteristic 






■V, 


fear 


vectors. Given this circumstance, the ct-sample principal components Y^ 
of are computed, as: ' 


c f = * <) (U [- M ? !; i = 1 > a > 3 - *• 


Then, the residual spectral signature of the specific pollutant is computed 


S“ C = u5 - [M? + 2 Y? e“ ]. (11) 

J J J i=l 1 ^ 

a? 55 a? p -,/p 

If [S| = [ 2 (S - ) ] x is significant relative to unity, then the pollutant 
3KL 0 

may he detected In sea water solely on the basis of its unique residual 

aC 

spectral signature. On the other hand, should [s| he negligibly small, then 

except possibly 'in a relative sense, a pollutant-dominated spectrum is 

J 

not uniquely distinguishable from "ordinary" ocean color spectra in sample a. 


EXPERIMENTAL RESULTS 


Upvelled radiance spectra IT. were measured with an Ebert spectrometer 

J 

aimed at nadir through the bottom of an NASA C-5& aircraft. The aircraft flew 
over the ocean at an altitude of 523 meters (lJOO ft), covering tracklines 
that included relatively clear offshore water, highly turbid Hudson River 
effluent, and pollutant dumps in the New York Bight Apex. The spectrometer’s 
field of view was approximately 30°, thus including a circle approximately 
270 meters in diameter. During the 1.5 seconds required to observe a single 
spectrum, aircraft motion displaced the field of view approximately 115 meters 
along the trackline. 


Hv I 
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A grating in an Ebert spectrometer reflects the diffracted spectrum 
through a slit. In this experiment, the slit was set to pass an approximately 
2.5 nanometer (nm) vide band of wavelengths onto a photomultiplier tube. The 
grating is rotated to scan the spectrum across the slit, so that a time record 
of the photomultiplier’s output voltage may be calibrated to yield radiance 
as a function of wavelength. For the MESA experiment, the GSFC Ebert 
spectrometer was calibrated using the diffuse integrating sphere at Goddard 
Space Flight Center. 

A selected subsample of the observed spectra was digitized into 55 

channels, each 5 mi wide over the wavelength interval 420 to 695 nm. This 

particular representation was selected to facilitate comparison with earlier 
12 

results * . The individual radiance spectra were then calibrated and 
normalized to obtain ocean color spectra as defined by equation (l). 

Five U. spectra measured in the vicinity of the acid waste dump site were 

d 

set aside for separate analysis. The remaining sample of 32 spectra 
(j = 1, 2, . , 55; n = 1, 2, . . .,32) was then subjected to 

characteristic vector analysis. 

The first four characteristic vectors (i — 1, 2, 3, 4; j - 1, 2, 

. .., 55) of the MESA sample covariance matrix are shown, together with the 
sample mean spectrum Mj (j “ 1, 2, . . . , 55 ) 3 as dashed curves in figure 1. 
Also shown in figure 1 are similar results for Mueller' Oregon data (OSU 
sample) and Hovis’s^ ocean color spectrum measurements made on his 1972 NASA 
Convair 990 flights (CV 990 sample). Inspection of table 1 shows that the 
first four characteristic vectors account for 99 percent of total sample vari- 
ance in each of the three samples. Internally, therefore, each of the three 


camples has exactly five 99 percent significant degrees of freedom associated 
with N and the first four principal components (Y^, Yg, Y^ 3 Yj,) . 

That the samples merely have the same number of 99 percent significant 
degrees of freedom is not, however, sufficient evidence to support the hypoth- 
esis that standard characteristic vectors are applicable to all three samples. 
Obviously the parameter IT is common to all three samples. The key question, 
then, concerns the degree to which the characteristic vectors are 

mutually interchangeable between samples. Intersample agreement can be 
measured by the orthogonal discrepancy terms and ^13 

equations (4) through (8), and by the percent of variance retention in 

P 

cross-sample representations (eq. (9) )* Pending results of a larger sample 
analysis now in progress, part of which will correct for intersample 
spectrometer calibration discrepancies, will be neglected and it will be 

assumed that 


6 iS 


le a I 

1 E i3 l 


k 

= 1 - { 2 
k=l 


«W 2} 


1/2 


( 12 ) 


Characteristic vector coupling coefficients C ^ p (eq. (5) ) are given 

in table 2A for a = MESA and 3 = OSU samples, and in table 2B for a = MESA 

and B = CY 990 samples. The magnitude of orthogonal discrepancies l E ^| 

as defined by equation (ll) a are given in table 3 S together with cross-sample 

variance retention V?. 

p 

ct 

As indicated by the large values of the diagonal elements (table 2) 

the MESA characteristic vectors are reasonably well alined with their 
counterparts for the OSU and CV 990 samples. The percentage variance 
distributions are also similar, although the proportion (about 90 percent) of 
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sample variance associated -with is larger than the corresponding propor- 

tion (approximately 80 percent) for the other tvo samples. This may "be due, in 
part, to the relatively restricted subsample on which the present MESA results 
are "based. 

Of more importance than either of the variance direction similarities 

shown "by tables 1 and 2, however, are the very small values of je^j in 

table 3. These values show that only a few hundredths of the total lengths of 

and are not completely transferable between any of tbe samples. And 

only about a tenth of and is not similarly transferable. Because 

of this overall cross-sample equivalence in the first four characteristic 

vectors, 96 percent to 97 percent of total sample variance is retained if MESA 

spectra are parametrised with OSU or CV 990 vectors, and vice versa. This 

result strongly implies that four standard characteristic vectors could well 

provide a common five percent significant degrees of freedom parametrisation 

basis for the three samples of ocean color spectra. 

Turning finally to the acid waste signature, two spectra from within 

r 

a recent acid waste dump and on UJ spectrum from 200 meters beyond the 

r 

visible edge of tbe dump were selected for analysis. The two acid U? spectra 

J 

r 

are shown as solid lines in figures 2A and 2B. The third in spectrum is in 

figure 2C. These spectra, together with the MESA mean and characteristic 

vectors, (fig. 1), were substituted into equations (10 ) and (ll) to obtain 

residual spectral signatures S/, which are shown as dashed curves in figure 2. 

1 

Obviously, the two acid U. spectra have a significant residual spectral 

3 

Ct£* 

signature S. , whereas the third spectrum does not. 

Visual comparison of the overall and residual acid spectra (figs. 2A and 
2B) with the reflectrance spectrum of iron hydroxide powder- (not shown) suggests 


that a precipitate of that substance, rather than of iron oxide, nay be 
responsible for the observed spectral signature of the acid dump. 

CONCLUSIONS 

The results of an ocean color spectrum analysis have been presented 
which use the method of characteristic vector. The following conclusions 
are made from the analysis: 

1. There is strong evidence to support the working hypothesis that a 
standard set of a four characteristic basis vectors plus a standard 
"mean" vector may be used globally to parametrize "ordinary 11 ocean color 
spectra. An individual spectrum would be parametrized by its mean 
radiance N and four standard principal components (Y^, Y^,, Y^, 
representing the five 97 percent significant degrees of freedom found in 
"ordinary" ocean color data. It Is emphasized, however, that much 
additional data is required before such a system could be accepted for 
standard use. By specifying "ordinary" ocean color it is meant to 
exclude waters containing industrial and acid waste materials. 

2. The acid waste material routinely dumped into the New York Bight 
appears to have a spectral signature with a unique component orthogonal to 
the characteristic vectors of ocean color spectra observed elsewhere. 

This unusual spectral signature may be caused by a FeOOH precipitate 

that forms after the pollutant enters the water. If constant from 
dump-to-dump and stable with time in the water, the residual spectral 
signature may be used to accurately discriminate the acid waste from, 
e.g., sediments or chlorophyll, and to man its dispersal patterns in 
remote ocean color imagery. 
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In closing, it is emphasized that the numerical results given here are 
preliminary and subject to modification in the future. The MESA and CV 990 
spectra were measured -with the same spectrometer, "but in different 
configurations and based on separate calibrations. The OSU spectra were 
measured with a different spectrometer which was independently calibrated 

by a different method entirely. Errors in any or all of the radiance calibra*- 

► ct ct 

tions of the spectrometers used will contribute to the error terms 

and 6. 0 . Recognizing these errors, a method of estimating the spectral 

ip 

inter calibration discrepancies between samples has been derived. 

Preliminary calculations suggest that using this method to adjust all spectra 
to a common calibration base virtually eliminates mean spectral orthogonal 
discrepancies . Corresponding improvements in and 6^ are 

expected. 
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Table 1. Percentage of total sample variance associated with 
the first four characteristic vector of sample convariance matrix for 
three independent samples of ocean color spectra. 



osu 

Sample 

C V 990 
Sample 

MESA 

Sample 

Id 

78 .6% 

82.9# 

90.85s 

! 2j 


15.1JS 

, 3.9% 

! 3J 

2 M 

1.2?» 

1 * 3 ^ 


0 , 9 % 

0.3% 

■ 0.7% 
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Table 2. Characteristic vector coupling coefficients C^g and , 
defining the projections of characteristic vectors e® . of sample a on 
those . of sample 3 and vice versa for MESA ocean coior data in comparison 
•with OSirand CV 990 samples. (See .text for explanation of terms). 

A. MESA {sample ct) vs OSU (sample 3), 


OSU (Sample 0) 


i k 

1 

2 

3 

1 + j 

1 

0.850 

-0.1+13 

O.Ol+O 

0.306 1 

MESA 2 

O.UiT 

0.861 

0,125 

i 

0.033 

(Sample a) _ 

D 

-0.170 

-0.11+7 

0.858 

0.198 j 

4 

-0.227 

0.168 

-0.132 

0.812 1 

B. MESA (Sample a) vs 

CV 990 (Sample 0). 


j 



CV 990 ( 

Sample 0) 

i 

. k 

1 

£ 

*n 

3 

1 + 

+ 

1 

0.933 

-0,255 

0.01+3 

0 . 21+1 ! 

1 

MESA 2 

0.288 

0.915 

0.044 

- 0.111 | 

(Sample a) 




j 

3 

-0.027 

-0.11+1 

0.836 

- 0.197 1 

■ i 

4 

-0.169 

0.199 

0.320 

0.746 



‘ •! 
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Table 3. Cross-sample unresolved orthogonality errors |£^d| aQ d- 
pereent of cross sample variance retention V9; for MESA ocean color spectra 
compared to QSU and CV 990 samples. (See text for explanation of terms.) 


Sample a 

osu 

MESA 

CV 990 

MESA 

Sample 3 

MESA 

OSU 

MESA 

CV 990 

KV 

0.006 

0.012 

0.002 

0.009 

i&i 

0.035 

0.020 

0.033 

0.019 


0.091 

0.122 

0.129 

0.102 


0.130 

0.109 

o.i4t 

0.184 


96.85? 

95. 1% 

9T M 

■ 97.135 
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Figure ] . - The mean vectors (M.) and first four characteristic . 

vectors (e..) of the sample covariance matrices for three samples 1 

of ocean color spectra measured in different geographical regions. - I 

Vectors M ■ are in units of relative radiance. Vectors | 

e . . are dimensionless and of unit length. I 

• J IS 

St jj 
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Figure 2.- Ocean color spectra (solid curves) and residual 
signatures (dashed curves) obtained as described in the main 
text in the vicinity of the acid waste dump in the New York 
Bight. 
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MULTICHANNEL OCEAN COLOR SENSOR (MOCS) 
by 

Gary V. Grew 
SUMMARY 

An algorithm for the remote detection and quantification of chlorophyll- A 
in coastal zone waters has been established for concentrations above 15 yg/l. 
Eor concentrations above and below this limit chlorophyll-A has been detected 
and identified by a signature extraction technique developed for analysis of 
data collected with (Multichannel Ocean Color Sensor), (MOCS) a multi spectral 
scanner. This paper presents MOCS data collected in the New York Bight in 
April 1975 as part of an NASA-NOAA Remote Sensing Program* Spectral signa- 
tures of algae, acid waste, and sediment are presented. 

INTRODUCTION 

It has long been recognized that the color of the ocean can be altered 
by the presence of algae. It is natural, therefore, to expect that man would 
investigate the possibility of remotely identifying and quantifying algae 
from aircraft and spacecraft. The desirability of generating contour plots 
of algae in the ocean is clear. Regions abundant with algae are of 
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particular interest to various user organisations for a number of reasons. 

For example, in order to support large populations of algae, regional -waters 
must be rich in nutrients. The sources and circulation of these nutrients 
are important directly or indirectly to the following groups: (l) to ocean- 

ographers in studying ocean and coastal currents (up-welling currents in 
particular), (2) to commercial fishermen in locating feeding grounds for fish, 
( 3 ) to ecologists in locating sources of pollution, and ( 4 ) to marine 
biologists in studying red tides. In the latter case, monitoring of red 
tides, dense populations of certain dinoflagellate algae, are important 
because they can he toxic to sea life and man and have been detrimental to 
the seafood industry, especially in. the states of Florida, Massachusetts , 
and California. Thus, interest in large area mapping of algae in the sea by 
remote-sensing techniques has been brewing since the onset of the U.S. 
space program. 

The feasibility of remotely detecting algae has been demonstrated by a 
1-3 

number of researchers. However, it appears that, as yet, no generally 
accepted algorithms for the identification and quantification of algae 
remotely has been established. The major difficulty is that the total 
upvelling light is dependent on so many variables , as illustrated in figure 1 
The variables shown in the figure are dependent upon such factors as solar 
elevation, atmospheric conditions, sea surface roughness, and composition 
and vertical distributions of the suspensoids in the sea, water depth, and 
ocean floor composition. Thus, an algorithm for algae that appears to 
work for one set of data may not work for a different set. An algorithm must 
be found that negates or minimizes the influence of all other variables. 


Fortunately, algae possess characteristic absorption spectra ■which can 
be used to some, as yet 01111010 ™ , degree. The ability of plant cells to absorb 
light depends upon their pigments. Figure 2 shows examples of absorption 
spectra measured by Yentsch.^ These curves show a strong absorption pealt 
at about 6 75 nm in the red region of the spectrum due primarily to chlorophyll - 
A and a second broader peak at ^50 nm in the blue region resulting from the 
combined absorption of chlorophyll and carotenoid pigments. The use of these 
two major absorption characteristics of plants is currently under investigation 
by the author for remotely quantifying alga.e in the sea. As will he discussed 
in later paragraphs, the chlorophyll-a peak in the red is useful at high 
concentrations (>15 yg/1 of chlorophyll-a) , but not at lower concentrations 
where the strong absorption characteristics of water (fig* 3) overshadow 
that of algae. At lower concentrations the absorption maxima in the blue is 
still usable, but it becomes increasingly more difficult to separate algae 
from sediment, the other major suspended matter in the ocean. It is the pur- 
pose of this paper to report the progress to date made in resolving these 
problems . 

The data presented in this paper were collected using Multichannel 
Ocean Color Sensor (MOCS), a multispectral scanner. This instrument was developed 
under NASA contract by TRW, Inc., as part of the Advanced Applications Flight 
Experiments Program (AAFE) * MOCS is being flown on an NASA Wallops C-54 
aircraft with the prime objective of collecting ocean color data that can be 
used in a computer program developed at Langley for extracting spectral 
signatures and algorithms that uniquely identify algae, sediment, and 
pollutants in our coastal waters. Results presented in this report were 
provided by MOCS data collected on an NASA-NOAA flight in the New York Bight 


on April 13, 1975. A good data set of chlorophyll spectra was collected 
on that mission. This data set is supported hy some ground truth information 
provided hy ITOAA. 

MOCS INSTRUMENT DESCRIPTION 

MOCS Is a visible imaging spectroradiometer ■which performs multispectral 
scanning electronically. It has no moving parts. MOCS was specifically 
designed for measurements of small differences in ocean color from space. It 
measures the intensity In 20 spectral hands at each of 150 spatial sites of 
the ocean across the field of view. MOCS is unique In that It uses only one 
detector and, as a result, it is compact and very light, weighing only 23 
pounds. 

Figure 4 is a schematic of the optical arrangement of MOCS and a listing 
of its specifications. The results of recent measurements of the center 
frequency of each hand ere given in table 1. In operation, light from the 
water is focused hy the objective lens on the entrance slit. The instrument 
is designed to form a high-quality optical image of the ocean surface on the 
slit, so that light from one edge of the field of view is imaged at one end 
of the slit, light from the center of the field Is imaged at the center of the 
slit, and so forth. The light is then collimated, dispersed by a biased 
transmission diffraction grating, and reimaged on the face of the image 
dissector. The resulting image consists of a large number of adjacent spectra 
each one composed of radiation coming from a different site across the 
instantaneous field of view. The spectra are scanned in sequence in a raster 
pattern on the photosensitive surface of the tube. The resulting video 
signal is a measure of the spectral intensities of the light coming from each 


of 150 spatial sites* The scan rate is such that the whole raster (one 
frame) is read out in the time (286 msec) it takes the spacecraft or aircraft 
to move forward over the ocean one resolution element. The scan is then 
repeated at a rate of about 3-5 frames/sec to give contiguous coverage of the 
ocean. MOCS has a 17*1 degree field of view and a spatial resolution of 
2 by 4 milliradians . 

The present MOCS system has an alternate mode of operation made possible 
by changing the image dissector. In this mode the spectral resolution (5 nm) 
is increased by a factor of 3 at the expense of a factor of 3 reduction in 
spatial resolution* 

The output of MOCS is fed to an A/D converter and stored on magnetic 
tape. The bit rate from the converter is 137 kbits/seconds. A detailed 
description of MOCS and its associated electronics can be found in reference 6 . 
Figure 5 is a photograph of MOCS. 

NEW YORK EIGHT MISSION 

Bemote-sensing experiments were conducted in the New York Bight between 

April 7-17 » 1975 j as part of cooperative program between NASA and the 

National Oceanic and Atmospheric Administration (NOAA). The objective of 

this mission was to evaluate multispectral, infrared, microwave, and 

photographic remote-sensing techniques for determining surface circulation 

features and for detecting sediment, chlorophyll, and other water quality 

7 

parameters in the New York Bight. This mission was one of several Marine 
Ecosystems Analysis (MESA) programs that are being planned by NOAA. The MESA 
program is focused on providing information to better understand and minimise 
man's impact on the coastal zone. 
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Aircraft flights of the MOCS Experiment -were conducted on April ICth 
and April 13th, She aircraft flight lines, environmental conditions, end 
vater masses were generally the same on both days, The surface winds and, 
therefore, the waves were higher on April 13th. Good MOCS data were collected 
on both days and analysis of the data revealed similar results. Since more 
extensive ground truth data were collected by HGAA on April 13th, that mission 
is presented here. MOCS and other experiments were flown on a Wallops C-54 
aircraft at 5*33 km altitude. The flight lines consisted of five parallel 
tracks about 3*35 km apart and 55 km long in northwest-southeast direction 
as shown in figure 6. These flight lines are superimposed on an ITOAA generated 
map showing (l) the locations of water samples collected by helicopter, 

(2) the chlorophyll-a concentrations measured from these samples, and 

(3) possible contour lines of the apparent algae bloom on April 13th. 

Three sampling stations near Sandy Hook indicated high concentrations 
of chlorophyll-a (fig. 6). Based on discussions with UOAA of their 
preliminary analysis of the ground truth data, it appears highly probable 
that the primary constituent in this plume was algae. Figure 7 is a plot of 
averaged MOCS data for spectral band 5 (460 nm) along line 1) beginning just 
offshore of Staten Island. A plot of the ratio of bend l4 (601 ran) and band 5 
for the same data is shown in figure 8 along with the approximate location of 

ground truth points collected by means of helicopter and the IIOAA ship 
Kelez. The chlorophyll-a levels and the times between each vatei sample 
pickup and the C-54 overflight are indicated. 

The ratio plotted in figure 8 is used by the author as an ’’indicator'* of 
turbidity. Selection of these bands was based on the examination of a large 
number of MOCS spectra which revealed that the maxima ar.d minima of the 


spectral variations of light from the sea were at or were near hands l4 and 5, 
respectively. Mote that the inflection point of the attenuation length 
curve in figure 3 occurs at hand ih. Under ideal conditions (clear, 'windless 
day and solar elevation 20° - 50°) this turbidity indicator (Tl) ratio using 
uncorrected MOCS data is about 0.28 for clear vater and 0.60 for very turbid 
vater . 

Figure 8 indicates turbid water in the lower bay between Staten Island 
and Sandy Hook, The plume boundary beyond. Sandy Hook is quite evident. 

Beyond the plume there is a relatively clear water region followed by a 
sulfuric acid waste dump. By inspection the correlation between the TI 
ratio and the chlorophyll- a levels is reasonably good. Figures 9 and 10 are 
plots of the TI ratio for all five flight lines on April 13th and April 10th. 
The turbid plumes and acid waste dumps are evident in both figures. 

In the following paragraphs a signature extraction technique is outlined 
and at the same time spectral signatures of the plume water and acid waste 
in the Bight are presented. 

SPECTRAL SIGNATURE EXTRACTION 

As previously stated, many environmental parameters influence the spectra 
of upwelling light from the sea. The contribution of the backscattered light 
from a particular suspensoid is often a small percentage of the total signal. 
Thus, spectral data collected over similar water masses but under different 
environmental conditions yield different spectra. It is } therefore, extremely 
difficult to extract any information from a single spectrum. An investigator 
must examine the variations in the spectral data that result from variations 
in the concentrations of a suspensoid. 
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The MOCS data used here to demonstrate the MQCS signature extraction 
technique is a set of spectra collected on April 13th along line D. This 
data set was collected in about 20 seconds over a 1 -mile path in a region 
where the chlorophyll-a concentration is assumed to have varied between 6 
and l4 jig/1 (see fig. 8 ), This assumption is based on the ground truth and 
MQCS data consistency. Figure 11 shows rather uninspiring sample reflectance 
spectra from this chlorophyll" data set. One needs to examine the spectral 
differences between these spectra. Figure 12 shows plots of the fractional 
differences between spectra selected at random from the data set. Although 
Interesting features appear, they are partially lost In the '’noise." A 
technique is needed to extract an average fractional difference spectrum. 

This can be derived in the following manner. 

In figure 12 there appear to be fair correlations between spectral bands. 
For bands 12 (568 nm) and l4 (601 nm) , for example, a linear regression 
equation can be determined for the chlorophyll data set of the form 


'12 


- A 


I lU + B 


(1) 


where I ? (I is the signal of band 12 (band l4) for a given spectrum. A 
and B are constants. For the data set, the correlation coefficient for equation 
(1) was found to be 0.94. Likewise, an equation and correlation coefficient 
can be determined for each of the other 18 MOCS spectral bands with respect 
to band l4* These equations can be written 


I. 

l 


“ A i,l^ X lU + B i,l4 


1 = 1,2 

i ^ l4 


. . 20 


( 2 ) 
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Nov the relative fractional change s g(i,l4),fcr band i for the curves in 
figure 12 can be wr'tten 

i dI - 

g(i s llO = 7--^-- (3) 

I. dI l4 

i 

or from equation (2) 


g(i s l4) 



i = 1, 2 . . .20 
i ^ 14 


(4) 


where 31 is the average signal of band i. Plots of g(i s l4) (normalized 
by the peak value of g(i s l4)) and their corresponding correlation coefficients 
for the chlorophyll data set are shown in figures 13 and l4, respectively. 
Compare figures 12 and 13- Figure 13 can be considered to be a variant 
spectral signature of the chlorophyll data set. Figure 15 shows plots of 
variant spectral signatures determined from 50 different subsets of 75 spectra 
each from the chlorophyll data set. The plot in figure 13 is actually an 
average g(i,l4) of the data in figure 15- 

Selection of band l4 in Fq. (4) was somewhat arbitrary. For a selected 
band j equation (4) can be rewritten in general as 

A. . 

g(i s j) = 

I. 


i = 1. 


i # j 


(5) 


A normalized variant spectral signature vising j — 12 is plotted in figure l6. 
Figures 13 and l6 are essentially the same. Note that the undefined 


quantity g(l2,12) has "been art ifie ally added to figure l6 by interpolation 
using the values of hands 11, 12, and 13 in figure 13. The variable Jk in 
equation (5) re-introduces the atmospheric scattering component of the upwelling 
light into the spectral signatures. This component , however, does not 
seriously influence the spectral features; it merely decreases the magnitudes 
of g(i,j), particularly in the blue bands. The term is used for two basic 
reasons. First, the algorithms selected for each suspensoid will probably 
consist of ratios of two or more spectral bands. Variations in these ratios 
depend on the fractional changes, Al^/l^ a of the total signal for each band. 
Second, g{i,d) is independent of the sensor used, that is, for all ocean color 
systems at the same altitude, the spectral signature of a particular water 
mass obtained by means of equation (5) would be identical. Ho sensor calibration 
is required. 

Further examples of variant spectral signatures exti’acted for the April 
13th mission are shown in figure 17- The data set for this case was collected 
over a sulfuric acid waste dump in the New York Bight. Differences between 
the chlorphyll and acid waste signatures can be seen in figure 18. 

ALGAE 

Additional algae signatures shown in figure 19 were extracted from the New 
York Bight data. One was extracted from a data set in which the chlorophyll- 
a concentration varied between l4 arid 2k pg/l along line D on April 13th 
(see fig.. 8). The second signature was obtained from the plume on April 10th. 
Other similar spectra have been extracted from data collected with KOCS in 
the Chesapeake Bay region. 
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A prominent feature of these signatures is the chlorophyll-a absorption 
peak at hand 19 (678 run) in the red region of the spectrum. Note that the 
absorption peak is negative going in the reflectance spectra (fig. 19) ^ 
positive going in the absorption spectra in figure 2, A possible secondary 
absorption peak due to chloroph^, 1- a or some other pigment is evident at 
band 1 6 (631 run). Absorption in the blue region of the spectrum due to 
chlorophyll-a and carotenoid pigments is indicated by the sharp slopes of 
the spectra beginning at band l4 (601 nm) and leveling off near band 9 
(521 nm). 

Since the relative signal of each band in figure 19 is positive, some of 
the variable component of the upwelling light is being backscattered rather 
than absorbed by the algae pigments. Whether this backscattering is due 
primarily to algae itself or some other process in the sea is unclear. 

Further distortion of the absorption spectra due to increased backscattering 
is evident in figure 20. This signature was extracted from a data set which 
included spectra all along line D from Staten Island to the plume boundary, 
excluding shallow water regions. A more extensive knowledge of the composi- 
tion of the suspended matter in the water is needed before these spectral 
signatures can be understood. 

Even so, the chlorophyll-a absorption in the red region of the spectrum 

should be useful to some extent for quant if ying algae. As the concentration 

of algae increases this chlorophyll-a feature becomes more pronounced. This 

principle is illustrated in figure 21 in which the ratio I /I versus the 

j- ( 

ratio of ^20^19 f° r points along line D. As the 

chlorophyll-a feature becomes more prominent with increasing concentration 
I19/I17 <3 - ecreases and ^po^l? increases as verified by the ground truth data. 


With, the exception of the "tail 11 in the lower right of the plot, a good 
linear relationship exists between these ratios, and should, therefore,' be 
useful in quantifying chlorophyll~a in the sea. Analysis of the data points 
in the "tail" indicates that this deviation is due to bottom reflection in 
very shallow water (<3 meters). Once the effect of shallow water is clear, an 
optimum algorithm for chlorophyll in very shallow water can be selected using 
spectral bands 17, 19s and 20. Perhaps a combination of the two ratios, 
I^p/Il^ and I 2Q /I 199 should be used. For example I 2Q • is plotted 

in figure 22 versus far comparison. Notice the effect of the 

shallow water data for each ratio . 

Since only a very few data points occur over very shallow water, the 
ratio was used in figure 23 to generate the false color map 

of ehlorophyll-a concentration along line D. The algorithm used for this map 
was 


IT = C 



( 6 ) 


where N is the chlorophyll-a concentration and C and P are constants. 
Notice the ground truth measurement of 32.2 yg/l in figure 6 at a station off 
Sandy Hook, but outside the field of view of MOCS for line P. It is in 
agreement with the color coded level, >_29 yg/l, indicated near Sandy Hook in 
figure 23, which was determined from equation (6). Figure 2k shows plots of 
^° r ^ ve ^-^Sh^ 011 April 13th. The algae plume appears to 

be quite strong in New York’s Lower Bay and along the New Jersey shore in 
agreement with the ground truth data in figure 6. Line A shows high concentra- 
tions of algae in Rockaway Inlet and much lower concentrations offshore of 
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Rockaway Beach as verified also in figure 25- A similar plume existed on 
April 10th as shown in figure 2 6 except that it extended farther seaward and 
its boundary was more sharply defined, , 

As demonstrated by the data in figures 21 and 25, the chlorophyll-a peak, 
at band 19 does not appear to be useful for quantifying algae much "below 
15 lig/1- As the algae population thins out, absorption of light by the ocean 
increases. A point is reached in which absorption of light by water is 
much greater than by the chlorophyll-a pigment at band 19 and the feature 
becomes lost in the noise. Thus, for low concentrations, another algorithm 
must be found - one using spectral bands in the blue and green regions where 
absorption by water is not as strong. Several promising algorithms are 
currently under investigation. 

SEDIMENT 

Two authors demonstrate elsewhere in this document that there is a high 
correlation between the sea truth measurements in the Bight of suspended 
solids and chlorophyll-a. This result practically precludes the use of this 
sea truth data to separate out these two suspended materials. Moreover, this 
high correlation suggests that the primary constituent in the plume was algae, 
as previously indicated. 

A variant spectral signature of sediment, compered in figure 27 with an algae 
signature, has been extracted from MOGS data on another mission (ref. 8). 

This sediment signature exhibits no unique features, suggesting that clear 
separation of sediment from other constituents in the ocean may be limited. 
However, two separation techniques are currently under investigation by the 
author. Their preliminary results are beyond the scope of this report. 


CONCLUDING- REMARKS 

Variant spectral signatures of algae, sediment, and acid waste have 
been obtained by a signature extraction technique developed specifically 
for analysis of MOCS data. The signatures were normalized such that they 
are independent of the ocean color sensor used; no instrument calibration is 
required. As demonstrated in this paper, a distinct advantage of the variant 
spectral signatures is that they can be used directly to establish algorithms. 

The chlorophyll-a absorption peak in the red region of the spectrum is 
clearly evident in the spectral signatures of algae. This feature appears 
to be useful for quantifying chlorophyll-a above about 15 yg/1 by use of bands 
IT (6^7 2a): 19 (678 rna), and 20 ( 69 k ran). A comparison between the ratios 
and ^20^19 clearl y separates out chlorophyll-a for the algae in 
the New York Bight on April 10th and 13th. The latter ratio can be used in 
deep water. For very shallow water (<3 meters) perhaps the ratio 

g 

IgO * I^/Cl.^) should be used. Whether or not these algorithms can be 
used for most saltwater and freshwater species is, of course, uncertain. 
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Table 1 MOCS Spectral Bands 
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Figure 1.- Interactions of light with the atmosphere and the ocean 
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Figure 2.- Absorption spectra of extracts of plankton algae 
(from Yentsch^). 
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Figure 4.- Optical schematic and specifications of M0CS 




Figure 5.- Photograph of MOCS 
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Fiqure 6.- Approximate flight lines of the NASA Wallops C-54 on 
April 13, 1975, superimposed on a NOAA generated map 
showing locations of stations where chlorophyl 1 -A 
(pq/l) measurements were determined from water samples 
collected by helicopter. 





FRAME 

Figure 8*- Averaged data points along line D for ratio of band 14 
(601 nm) and band 5 (460 nm). Chlorophyll-a levels 
(ug/1) collected by helicopter and the NOAA ship Kelez 
are indicated along with times (hours) in parentheses 
between water sample pickups and C-54 overflight. 
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3*- TI ratio for all five flight lines in the New York Bight 
on April 13, 1975. 
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Figure 10.- TI ratio for all five flight lines in the New York Bight 
on April 10, 1975. 
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Figure 11. -Sample reflectance spectra from MOCS data collected at 
5.33 km altitude along line D. 
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Figure 15.- Variant spectral signatures of 50 subsets of the chlorophyll-a 
data set, j =14. 
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Figure 18.- Comparison between variant spectral signatures of algae and acid 
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Figure 19.- Variant spectral signatures of algae from New York Bight data sets. 
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Figure 20.- Variant spectral signature of a data set consisting of spectra 
all along the plume water from its outer boundary to Staten 
Island. 
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Figure 23.- False color map of chlorophyll-a concentration along line D 
in the New York Bight. Arrows point to approximate locations 
of around truth stations. 
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REMOTE INFRARED MEASUREIMTS OF WATER THERMAL ANOMALIES DURING THE 

NEW YORK BIGHT EXPERIMENT 

Ly 

Charles A. Hardesty 
SUMMARY 

This paper presents thermal data collected in the New York Eight in 
April 1975 from the C-5^ aircraft platform as part of a KASA-NOAA cooperative 
program. These thermal data were collected for correlation with KOCE data 
and are not to he confused with thermal data collected onboard other flight 
platforms during this effort. Thermal data are presented here and compared 
with MOCS spectral signatures for algae. 
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INTRODUCTION 

* 

Remote measurement of sea surface temperatures using infrared radiometers 
has been widely studied and the technique used on many occasions. Absolute 
accuracies of better than 0,5°C are obtainable from aircraft platforms 
provided the infrared radiation measurements have been processed to remove 
the atmospheric effects and by using ground truth calibration. The infrared 
experiments flown on April 10th and April 13th were a part of the mission to 
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evaluate multispectral, infrared., microwave, and photographic remote sensing 
techniques for determining surface circulation and water quality features in 
the New York Bight. 

The purpose of including the infrared thermal scanning experiments on 
this mission (C~5t aircraft) was to correlate the observable thermal anomalies 
in the water with the remote measurements collected by the Multichannel 
Ocean Color Sensor (MOCS). Water is essentially opaque to the infrared 
radiation in the 8 to it micron range. Therefore, the upper* 0.1 mm of sea 
surface determines the radiance variation from point to point and. the 
temperature of the marine background. The thermal data collected have not 
been processed to remove the effects of atmospheric constituents. These 
data are presented in this report as temperature changes along the flight 
lines. The poor quality of the thermal scanner film records arid the magnetic 
data tapes prevented digitising the data and providing image enhancement 
ana false color maps as was originally planned. 

DESCRIPTION OF INFRARED INSTRUMENTATION 

A Precision Radiation Thermometer (PPxT-5), table 1, and a Reconofax IV, 
sable 2, were flown on the New York Bight mission. The FRT-5 consists of 
an optical unit and an electronic unit. The optical unit compares the amount 
of radiant energy emitted by the target with that emitted by an internally 
controlled, optically chopped temperature cavity, thus providing absolute 
values of infrared radiation levels. The unit uses a hiper immersed thermistor 
detector with the area of the active element only 50 micron square. An 
objective lens and special interference filter are used to limit the pass band 
to the atmospheric 8 to it micron window. With the 2 degree (2°) field of view, the 


instrument covers a ground swath of 15-9 meters and 186 meters for altitudes 
of 0.45 and 5-33 km (l400 and 1T S 500 ft.), respectively. 

The He conofax instrument shown in figure 1 is a single channel infrared 
mapper. Incoming radiation from ground targets is reflected from the *4 5 
degree surface of the rotating (11,000 rpm) scanner mirror into a parabolic 
mirror. The parabolic mirror focuses this energy via a small right angle 
folding mirror onto the sensitive surface of the detector. The varying 
intensity of the received energy produces an electrical signal. The 
resulting video signal is used to modulate a crater lamp and record the 
scanned image on TO mm film. A mercury doped germanium sensor cooled to 
approximately 30°K was used on this mission with an 8 to lU micron atmospheric 
filter. Since this Reconofax did not have an onboard calibration system, only 
relative infrared radiation levels were obtained with the scanner. 

The remote infrared and MOCS experiments were mounted on the Wallops 
Flight Center's C~54 aircraft platform and flown on April 10th and April 13tli. 
Flight data wex^e collected at the 0.45 and 5*33 km (l400 and IT , 500 ft) 
altitudes. The flight lines consisted, of five parallel tracks (A thru E.) 
about 3*55 km ( 1.92 nautical miles, n. mi.) apart and 55 km (2 9-9 n. mi.) long 
in the northeast-southeast direction a. vm in figure 2. Flight lines 1, 2, 
3 5 ^ , and the extended C--1 line, were f_ i on April 10th at an altitude of 
0.45 km (1400 ft. ). 

U0AA provided the ground truth data for both days, but the most extensive 
set of temperature data were collected on April 13th. In most cases there „was 
a considerable time spread between the NOAA ground truth data points and the 
UASA remote sensing data obtained, from the C-54 aircraft. 


DISCUSSION 


Figure 3 provides an excellent sample of the quality of the Eeconofax 
IV thermal scanner film records obtained during the spring mission. This 
figure shows flight lines A-l and D-l flown on April 10th at an altitude of 
5-33 km (17,500 ft,). Some of the fine detail has been lost in reproduction, 
but close observation of the film records shows several small ships (small 
dots) along with their cooler wake trails. This observed temperature 
difference is a result of cooler water being brought to the surface by the 
mixing motion of the moving ships. The thermal boundary off shore of Long 
Island is quite pronounced in this picture which indicates a temperature 
difference on the surface due to the mixing motion (flow) of the water. The 
sewage sludge dump is also quite apparent in the middle of the photograph and 
the presence of the wake from the dumping ship indicates that this was a fresh 
sewage dump. The sewage dump was first noted in flight line B-l and it can still 
be seen in run B-2, approximately 100 minutes later. The reason for this 
strong sewage sludge signature and its persistence has yet to be determined. 
Additional information is needed on the sludge temperature arid composition 
at the time of discharge to clarify this point. It is possible that the 
sludge enrichment of the overlying waters may have resulted in formation of 
extensive surface films. The old acid dump is not visible in the photo- 
graphs. This is reasonable, because the Eeconofax IV can only see a 

* 

temperature change of 0.3°K or larger in the top 0.1 mm surface of water, 
or a change in the emissivity of the surface. Other features to be noted or, 
the Eeconofax records are the small hairline scratches running parallel 
with the pictures along with two heavier dark lines. These marks were caused 
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by the film magazine in the scanner and makes it very difficult to digitize 
and analyze the data using a densitometer. 

Figure 4 is an example of the Reconofax records for the 0.45 km (1400 
ft.) low altitude radial lines 1 and 2 (reduced by 2/3 to permit printing on a 
single page) flora on April 10th. Line 1 starts Just South of Sandy Hook 
and ends at the bridge on Rockaway Beach. Shora are thermal gradients 
located approximately midway (arrows) between the land points, and just off 
the shore of Rockaway Beach. The PRT-5 low level record in figure 5 is on 
example of the last 15 minutes of data for the 85 km (46 nt mi.) "C" line 
on April 13th showing two thermal boundaries in the surface of the water. 

These boundaries correlate with the scanner film record and are verified 
by ground truth dara points 23 5 24, and 54. The Keles data points are 
approximately 1.4 C warmer than the airborne measurements which is reasonable, 
since the data have not been corrected for atmospheric effects as mentioned 
and the time spread between the data sets varied from 0.5 to 2.5 hours. 

The magnitude of the thermal anomalies at the boundaries varied between 
0.3 and 0.6°C. 

One interesting phenomenon was observed on the A-2 line flown at 5*33 km 
(17,500 ft) altitude on April 10th , and it appears on both the Reconofax 
(fig. 6 b) and FRT-5 (fig. 7b) records. This feature appears as many 
small thermal anomalies on the Reconofax: film record with a faint leading 
(boundary) edge to the anomaly area. Because the footprint (swath width) 
of the PRT-5 is 186 meters (6l4 ft.) for the high altitude pass, it tends to 
smooth or average the small gradients seen on the Reconofax record ( 16 - 
meter IFOV) for a total change of approximately 1°C over the area. The 
thermal feature mentioned was not apparent on the A-l (fig. 6 a and 7a) 


line flown approximately 30 minutes earlier, "but the gradient is still faintly 
visible in the B-2 line flown 30 minutes after the A-2 line. These high 
altitude flight lines were flown in the sequence of C-l, 1-1, D-l, A-l, C-2, 
A-2, D-2* and B-2 . The phenomenon first appeared on the A-2 line, hut the 
area was obscured by cloud coverage on the C-2 and B-2 flight lines , which 
prevented pin-pointing the time of the occurrence. This anomaly occurred at 
the approximate time of maximum current movement for this area of the Hew 
York Eight, thus indicating a possible convergence and mixing of the water 
mass. In some respects, the anomaly seen on the Reconofax record (fig. 6b) 
resembles the edge of the plume area indicated in figure 8, a sketch made 
from previous 1973 flight data over the Bight area. The MOCS record for the 
A~2 line also shows an increased amount of turbidity for this time frame, 
thus indicating in this case some correlation between the observed surface 
thermal anomalies and the MOCS data. 

Figures 9 and 10 are plots of the PRT-5 thermal data for the first five 
high altitude flight lines flown on April 10, and April 13, 19^5, over the 
Hew York Bight. These data can be compared to the MOCS data for retie 
in figures 24 ana 2b in reference 1. A comparison of these data sets 
suggests there may be some correlation between the surface thermal anomalies 
and MOCS data. 

C OK CLUE IMG REMARKS 


Thermal records for the FP.T-5 end Reconofax IV were obtained for all 
flight lines flown on. the C-5b aircraft for April 10, arid April 13, 197*3. 
The Reconoiax records for April 1 3th show less uucrmal activity than on 
April 20th* The quality of the Reconofax film records ard magnetic tapes 



prevented digitizing the data as originally planned. A comparison of the 

thermal IR and MOCS data sets indicated correlation in some cases, hut not in 

others. More consistent correlation between data sets may he obtained by 
+ < 

using different ratios of MOCS data other than the ratio. The thermal 

patterns seen on the Reeonofax records for this mission are indicative of 
water flow conditions and are similar in appearance (ref. 2) to the results 
obtained by the Environmental Research Institute of Michigan in their work 
over the Hew York Bight in 1973- 
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'■ TABLE 1, - PREOSffl RADTATTON THFRMOMFTFR (PRT-5) 

Sill A LIGHTWEIGHT/ PORTABLE BATTERY-POWERED THERMAL INFRARED RADIOMETER 
i SINGLE BANDPASS = 8 TO 14 «M WAVELENGTH 
i 2° FIELD OF VIEW OR 186 METERS a 5/335 METERS 
16 METERS a 457 M, ALTITUDE 


VO 

-p- 


PHYSTCA1 PARAMFTFR MFASIIRFD: THERMAL INFRARED RADIATION IN THE 8 TO 14 MICRON 


WAVELENGTH BAND TO ±0.5° C 

FARTH RESOURCES PROGRAM APPI .1 CATIONS; PROVIDES A TARGET TEMPERATURE REFERENCE 
MANUFACTURER; BARNES RADIATION 

PHYSTCAI SPECIFICATIONS; TWO MAJOR COMPONENTS - OPTICAL UNIT 

ELECTRONICS 


REQUIRES ON-BOARD RECORDER 
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TABLE 2, - 


SINGLE CHANNEL INFRARED LINE-SCANNING IMAGING SYSTEM 


THAT RECORDS RADIANT ENERGY IN THE 8 TO IN vM RANGE 


SPFCTRAI CHARACTERISTICS: 


SINGLE CHANNEL - 8 TO 14 nM WAVELENGTH 


■SPATIAI CHARACTERISTICS 


120° SCAN (±60° FROM NADIR) 3 MRAD INSTANTANEOUS FQV 


RESOLUTION! 


FOV a 5,335 M ALT, : 16 METERS 


FOV 8 457 M ALT.: 1.4 M 

GROUND SWATH a 5,335 Mj 18,300 METERS 

GROUND SWATH a 457 Mj 1,500 METERS 

PHYSICAL PARAMETERS MEASURED: INFRARED RADIANT ENERGY IN THE 8 TO 14 nM 

WAVELENGTH BAND TO +0,3 ° K 


IMAGE RECORDED ON 70 MM FILM 
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Figure 5 - PRT-5 RECORD - APRIL 13, 1975 - LOW LEVEL LINE C-l FLOWN AT AN ALTITUDE OF 457 METERS 
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MAPPING IKS HUDSON RIVER PLUME AND 
AN ACID WASTE PLUME BY REMOTE SENSING IN 
THE NEW YORK BIGHT APEX, APRIL 1975 

tv 

Herbert W. Johnson 
INTRODUCTION 

Suspended sediment and chlorophyll~ja concentrations and their distribu- 
tions have been identified as important environmental parameters for monitor- 
ing "water quality and pollution in coastal zones. Suspended sediment has 
been recognized as a natural tracer that may be used to measure flow and 
distributions in a water body, thereby providing information on pollutant 
concentrations and dispersions. Chlorophyll-a_ concentrations are a measure 
of the nutrient load and an indicator of current state of health of a water 
body. In addition, as yet unidentified materials in spectral anomalies such 
as plumes resulting from ocean dumping of acid wastes may be used to identify 
and map those features. Synoptic distributions of these environmental 
parameters may be determined from remotely sensed data, thereby providing 
information not readily available by any other means. Calibrated regression 
equations have been used by Johnson (ref. 1, 2) to quantitatively map dis- 
tributions of water quality parameters such as suspended sediment and 
chlorophyll-Ei from remotely sensed multispectral scanner data. 

It is the objective of this investigation to apply the methodology of 
reference 1 to a coastal zone environment and to develop calibrated regression 
equations to quantitatively map distributions of water quality parameters 
(e.g., suspended sediment end chlorophyll-a) measured in the sea truth 
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program.* In addition, other spectral anomalies in the scene, such as acid 
waste plumes, will he qualitatively mapped. 

EXPERIMENTAL METHOD 

Remotely sensed data were analyzed and interpreted using concurrently 
collected sea truth information. Data collections were made on April 9, 13, 
and l4, 1975 j with the most complete set collected on April 13. Results 
presented here will he limited to analysis of remotely sensed data collected 
on April 13, 1975 (ref. 3). 

Remotely sensed data were collected hy the Ocean Color Scanner (OCS) and 
Mitchell-Vinten cameras from a U-2 aircraft platform at a flight altitude 
of 19.7 kilometers (km) (65,000 ft.). The OCS is a 10-hand multispectral 
scanner with 10 hands in the visible and near IS spectral range. Band center 
wavelengths are from 433 nanometers (nm) to 772 nm for the nominally 20 nm 
wide hands , table I. Ground spatial resolution was about 75 meters (248 ft.) . 
Data were recorded onboard the aircraft on magnetic tapes with four hands 
(2, 4, 5 and 8) digitized during flight and the others recorded in an analog 
format. Initial data reduction and formatting were performed hy the NASA 
Goddard Space Flight Center (GSFC) with calibrated digitized data in Bands 1 
through 8 supplied to the Langley Research Center on computer compatible 
tapes (CCT). 

Sea truth measurements for this experiment were collected by helicopter 
over the 25 station matrix in the New York Bight apex. In addition, six 

*1 would like to acknowledge the contributions of Mr. Terry A. Nelsen of the 
AOML/NQAA for the collection and analysis of sea truth data during the 
experiment, and also, those of Messrs, Gilbert S. Balm and Robert M, Glasgow 
of the Vought Corporation who contributed to the computerised data analysis 
and mapping products. 
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samples were collected at IT sites of opportunity.” Sample locations are shown 
in figure 1. Stations 4 , 5 , 10 t 15 , and 21 were outside of the OCS scan and 
X4 and X5 could not he located due to missing LORAN data. Suspended sediment 
and chi or ophyll-ja determinations were made at most of the sample locations. 
Table II is a listing of suspended sediment and chlorophyll-ji concentrations 
analyzed in this report. In addition, particle size distributions were 
obtained within 24 hours hy Coulter Counter from water samples taken at these 
stations (fig. 2), Station locations for the sea truth sampling were 
determined by LORAN A with an estimated absolute accuracy of TO meters (220 ft). 

DATA ANALYSIS AND RESULTS 
Data Preprocessing 

Representative radiance values corresponding to the sea truth measurements 
were determined by locating the sampling station as nearly as possible, and 
then by taking the average of an 11 by 11 field centered at that location to 
obtain the representative value. This field size was empirically determined 
as the area required to compensate for uncontrollable spectral and spatial 
"noise 1 ’ or uncertainty. Radiance values (mw/cm^-ster-ym) associated with 
the sea truth sampling stations are listed in table III. 

Quantitative Data Analysis 

Stepwise Regression Analysis (SURA) was used to determine calibrated 
regression equations for quantitatively relating the sea truth measurements 
to remotely sensed data, as in references 1 and 2. In an SWRA, the inde- 
pendent variable (radiance in one OCS band) is selected that has the highest 
correlation with the dependent variable (water quality parameter; e.g., 
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suspended sediment or chlorophyll-^) * A number of variables are then 

selected consecutively until all of the independent variables (bands of OCS 

data) that make a significant contribution to determining the dependent 
* 

variable .are included in the 'regression equation and the others are ’'outside" 
the regression. Limiting the regression equation to significant variables 
reduces the anal ysis time and improves the accuracy of the results. The 
criterion for inclusion of variables is a 90-percent confidence level as 
determined by the statistical "F" test (see p. 171, ref. h for a discussion 
of the SWRA) . * 

Quantitative Analysis Results 

All values obtained in the sea truth program (that could be located in 
the scene) were included in the quantitative analyses of suspended sediment 
and chlorophyll-^ except those in the spectral anomaly associated with the 
acid waste dump plume (e.g„ all stations except 18, 23, Xk and XJ5 for suspended 
sediment and 18, 23, XI, X2, X3, Xt , X5 and X6 for chlorophyll- a ) . There 
were a total of 22 observations for suspended sediment and 18 for chlorophyll- 
a. The plume from the acid waste dumps will be discussed under qualitative 
mapping in a later section. 

Results of the SURA applied to suspended sediment and subsequently to 
chlorophyll-a. were as follows: 

Water OCS Band,RN Standard Correlation Correlation Range of 

Quality in Regression Error of Coefficient to Suspended Sea Truth 

Parameter Equation Estimate Sediment Measurements 

Suspended 

Sedi- 


ment E6 

1.39 

0.79 

- 

0 ,kS - 8,38 mg/£ 

Chloro- 

phyll- 

a R3» R6 

3.87 

0.83 

0.90 

2.20 - 24.30 mg/m 3 
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where M is the radiance in OCS band ft (i.e., E6 is radiance in band 6)5 
standard error of estimate is a measure of the scatter about the fitted 
regression line; correlation coefficient is a measure of the relative change 
among variables; correlation coefficient to suspended sediment is the linear 
correlation of that variable to suspended sediment ; and range of sea truth 
measurements are for the water quality parameter being analyzed. 

Comparison of the remotely sensed (calculated from the regression 
equation) and measured sea truth values for suspended sediment and chlorophyll 
^concentrations are shown in figures 3 and respectively. Deviations from 
the fitted regression line occur approximately randomly, thus the linear 
model appears adequate. Johnson (ref. 1* 2) previously indicated linear 
responses in this range for suspended sediment and chlorophyll- a.. 

Interpretations of the remotely sensed data should consider the high 
correlation between changes in suspended sediment concentrations and changes 
in chlorophyll-a concentrations. In this case s where the correlation coeffi- 
cient (0.90) of sea truth (suspended sediment to chlorophyll-a) is higher 
than for sea truth to remotely sensed data for either parameter, it should 
be interpreted as a combined response where neither parameter can be 
unamb igiously separated by the analysis. The high correlation between 
suspended sediment and chlorophyll -a_ is shown in their comparison in figure 5* 

Quantitative Mapping of Suspended Sediment and Chlorophyll- a 

Quantitative mapping of water quality parameter concentration distribu- 
tions may be determined from the regression equations. For each water quality 
parameter, concentrations are determined at each pixel (or equal spacings 
of lines or columns), this field of data is typically smoothed to remove 




?^,L‘, .r -tiyrcm. .i E-^vgiF Wi^L-iaMs* ■* 
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local spectral and spatial noise features and then a contour map is 
developed by a computerized plotting routine. The smoothing routine used 
in this analysis is an averaging on a line-by-line and eolumn-by- column 
basis in the data field where. the middle value is replaced hy the mean of it 
and the two adjacent values. Edge values remain the same. In this analysis 
each pixel in every third line of data was used to generate the field of data 
for mapping, Two smoothing passes were made for suspended sediment concen- 
trations and two were used for chlorophyll- a concentrations. 

Maps of suspended sediment and chlorophyll-a distributions, figures 6 
and 7.» also indicate i,he high correlation (0.9G S see fig. 5) between these 
parameters. The Hudson Biver plume extension into the Eight Apex s then a 
drift southward along the New Jersey shore s is clearly indicated in both 
figures. The previously located acid waste, plume is shown in both the 
suspended sediment and chlorophyll- distribution maps. Data interpretation 
of this feature is discussed in the. following section. The more concentrated 
portion of the acid waste plume is at station 18 which is in or near the 
legal dump area. Less concentrated areas of the acid waste plume are 
probably due to dispersion and to earlier dumps. Acid waste dump plumes 
have been observed to persist for greater than 2k hours (ref. 5). 

Qualitative Analysis of the Acid Waste Plume 

There are not sufficient sea truth points in the acid waste plume tc 
apply Quantitative analysis techniques* However s investigation of spectral 
changes associated with plume dispersions indicates that the spectral range 
most indicative of plume . changes are in the OCS band 3 to 6 spectral range. 
This is the same spectral range that indicates differences in suspended 
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sediment and chlorophyll- a. Thus, apparent changes in suspended sediment 
and chlorophyll- a. concentrations in the acid mate plume area should he 
interpreted as qualitative mapping of the acid plume (qualitative since 
materials in the acid plume are known to he different from those in the other 
locations of the experimental area). 

Spectral responses in the acid waste plume at locations of apparently 
different concentrations are shown in figure 8. For comparative purposes, 
spectral responses for the Hudson River Plume and background (Bight Apex) 
water are shown in figure 9* Radiance values in the acid plume are not as high 
or as spectrally wide as in the combined suspended sediment -chlorophyll-a. 
in the Hudson River plume. This is shown in the ratios of plume to ocean 
water radiances in figure 10. 

Particle Size Distributions 

Particle size distributions for the basic areas in the sea truth and 
remotely sensed data are shown in figure 2 for the Hudson River plume, acid 
waste plume and background water. In the Hudson River plume, the dominant 
size particles are small (6 phi(<£) or 15 micrometers (pm)); those in the 
acid plume are relatively large particles (3. 5$ or SO pm); and those in the 
background water are essentially bimodal (both large and small) in character. 

It should be emphasized that these are relative counts of particles and 
only relative distributions may be compared since little is known about 
the shape or specific gravity of the constituents. 

Identification of Pollutant Plumes 

Identification of the two major spectral features (e.g.,the Hudson 
River and acid waste plumes) would be possible in future experiments by 
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spectral responses in ‘the remotely sensed data. Figures 8 and 9 are 
radiances (mw/cm^-ster-ynO in the OCS hands for ranges of concentrations of 
the Hudson River end acid waste plumes, respectively. Note the general 
spectral similarity. However, when the highest concentration responses of 
the plumes are ratio ed to an in-scene calibration (ocean water) the 
differences are quite apparent above wavelengths of 580 nm (fig. 10 ). 

It appears that both radiance values and ratios to an in-scene source will 
be useful for pollutant identification and monitoring by remote sensing. 

Identification of plumes by sea truth measurements will be aided by 
the particle size characterization, figure 2. As indicated previously the 
Hudson Biver plume has predominantly small particles, the acid waste plume 
has predominately large particles and the ocean water is bimodal in character. 

CONCLUDING REMARKS 

Remote sensing may be . effectively applied to locate, identify and map 
major water features in the New York Bight Apex, A multispectral scanner 
from an aircraft platform was used to map the synoptic dispersion patterns 
of plumes from the Hudson River and ocean dumping. The ocean dumping plume 
was from acid waste dumps in the Apex on April 13, 1975. . Results indicate 
that the types of plumes may be identified from remotely sensed spectral 
responses without concurrent sea truth. In addition, particle size distri- 
butions from the sea truth measurements provided information to further 
characterize the river and acid veste dump plumes. 

Digital data analysis techniques were used to develop quantitative 
relationships between remotely sensed data and the two water quality 
parameters of interest, suspended sediment and chlorophyll- a_. Due to the 


high correlation in the measured values of these two parameters in the* sea 
truth program 9 it was not possible to determine unambiguous distributions 
of these parameters; however, their mutual changes have heen used to map the 
dispersion of the Hudson River plume. 

Qualitative and quantitative mapping of pollution related features were 
made for the Hew York Bight Apex. Additional sets of data with expanded 
sea truth measurements will he required to further evaluate the analysis 
techniques applied to this set of remotely sensed data. 
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TABUS I SPECTRAL CHARACTERISTICS OP THE OCEAN COLOR S CAREER 


TEE CHANNEL SCANNER WITH 75-METER RESOLUTION, 
SPECTRAL BAND: 

433 NANOMETERS + 10 run, 

* 471 NANOMETERS + 10 m, 

509 NANOMETERS + 10 mi s 

* 547 NANOMETERS + 10 nm, 

* 583 NANOMETERS +10 mi, 

620 NANOMETERS + 10 nm, 

* 662 NANOMETERS + 10 m, 

698 NANOMETERS £ 10 nm, 

733 NANOMETERS + 10 nm, 

772 NANOMETERS + 10 nm, 

* DIGITIZED ONBOARD U2 S 
INTERNAL CALIBRATION ONCE PER SCAN LINE. 
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TABLE II SEA TRUTH MEASUREMENTS* 


STATION 


1 

2 

3 

6 

T 

8 

9 

11 

12 
13 
ll 
1 6 
IT 

18** 

19 

20 
22 

23 ** 

2 b 

25 

XI 

a£ 

X3 

x6 


SUSPENDER SEDIMENT 
CONCENTRATION 
nig/1 

6. lib 
3.210 
1.700 

а. 380 
2.710 
2.310 

1.100 

б. 290 
2.550 
1.930 
1.060 
3.600 
0.1*60 
9.030 
0.680 
O.56O 
0.710 
1.69C 
1.360 
0.690 
1.730 
0.790 
0.580 

5.100 


*0nly date., in tie 0C3 scar, are shown. 


CHLOROPHYLL A 
COK CM4TRATIGK 

irg‘/m- J 

15.0C0 

1.600 

5.300 

21.300 

6 . 100 
6.600 
5.000 

15o00 
1 1.000 
1.600 

3.000 

17.800 

2.700 

2 . lOO 
2.2C0 
1.2CG 

2.100 
5.300 
2 . TOO 


**In acid waste plume* 
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TABLE III REPRESENTATIVE RADIANCE VALUES AT THE SEA TRUTH STATIONS 


Station OCS BAND RADIANCE 

mw/ cm "ster-pm 

Band 



l 

2 

3 

4 

5 

6 

7 

8 

1 

23.94 

19.68 

14.52 

11.72 

9.16 

6.66 

5.00 

3.00 

2 

24.57 

20.35 

14.89 

12.09 

9.36 

7.01 

5.26 

3.37 

3 

23. 74 

19.93 

14.68 

12.07 

9.61 

7.04 

5.08 

3.09 

6 

23.77 

19. 66 

14.51 

12. o 4 

10.49 

7.98 

5*90 

4.00 

7 

23.54 

,19.62 

l 4 . 4 i 

11.98 

9.82 

7.00 

5. 00 

3. 00 

8 

23.08 

18.67 

13.45 

11.00 

8.88 

6.00 

4.17 

3.00 

9 

23-54 

19.70 

14.53 

11.08 

9.00 

6.19 

5. 00 

3.00 

11 

23. 14 

18.89 

13.58 

11 . 46 

10.00 

7.00 

5.00 

3.10 

12 

f 

23.46 

19.36 

13.94 

11.28 

9.01 

6.60 

5.00 

3.00 

13 

23.42 

19.45 

14.00 

11.12 

9.00 

6.12 

5.00 

3.00 

Ik 

22.90 

18.37 

13. 40 

10.38 

8.00 

6.00 

4.00 

2.94 

1 6 

22.96 

19.34 

14.25 

11.45 

9.88 

7.07 

5.07 

3.30 

IT 

22.44 

18.08 

13.32 

10.72 

8.17 

5.98 

4.00 

2.98 

19 

22.83 

18 . 49 

13.47 

10.93 

8.12 

6.00 

4.17 

3.00 

20 

22.75 

18.92 

13.45 

10.08 

8.00 

6.00 

4.00 

3.00 

22 

22.31 

18.05 

13.51 

lo. 4 i 

8.12 

6.00 

4.06 

3.00 

24 

23.05 

18.66 

13.45 

10.47 

8.02 

6.01 

5.01 

3.00 

25 

22.71 

17.93 

12.46 

9.50 

7.03 

5.41 

4.00 

2.81 

XI 

23.05 

18.79 

13.59 

11.15 

8.58 

6.02 

4.84 

3.00 

X2 

22.12 

18.15 

13.93 

11.81 

9.23 

6.00 

4.61 

3.00 

X 3 

22.19 

18.12 

13.86 

11.51 

8.89 

6.00 

4.47 

3.00 

x6 

23.73 

19.68 

i 4 . 4 t 

ll. 4 l 

9.04 

6.38 

5-00 

3.00 
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Figu re 1. - Location of sea truth stations in the New York Biciht Apex 
on April 13, 1975, 
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Figure 2. - Particle size distributions for the Hudson River and 
acid waste plumes, and N.Y. Bight Apex ocean water 
on April 13, 1975, 
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Figure 3, - Comparison of measured and remotely sensed values 
of suspended sediment 
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Figured - Comparison of remotely sensed and measured values 
of chlorophyll a . 
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Figure 5. - Comparison of measured values of Chiorophiyy a and 
suspended sediment at sea truth stations. 
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Figure?. - Quantitative distribution of chlorophyll a (mg/m 3 ) in N. Y. Bight 
Apex April 13, 1975. 
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Figure 8. - Spectral responses in the plume resulting from acid 
waste dumps in the New York Bight 
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Figure 9. - Spectral responses of Hudson River plume at different 
concentrations and Ocean water on April 13,1975. 
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Figure 10. - Ratios of pollutant plume to ocean water radiances 
in the New York Bight on April 13, 1975. 
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SUMMARY 

Comparison o? broadband black and white photography at near-infrared 
wavelengths with green-yellow wavelengths shows that a substance being 
dumped in the Rev York Bight sludge dump area has a different spectral 
signature relative to the water background than does a substance from the 
acid dump area* Relative to the background water, the substance being 
dumped in the sludge dump area has greater reflectance at near-infrared 
wavelengths than at green-yellow wavelengths. The substance in the acid 
dump area has greater reflectance than the background water through both 
optical filters. In addition, the reflectance from the substance in the 
acid dump area is much greater in the green-yellow spectral region than it 
is in the near-infrared region. 

Film density data, obtained from the green-yellow photography, show that 
waste dumping in the Hew York Bight increases the backscattered sunlight 
from the water over a much greater area than the designated dump site areas. 
It is clear from the photography that acid waste is the primary substance 
responsible for the increased backscattered sunlight, but the green-yellow 
wavelength data suggests that sludge dumping is also contributing. 

Film density data indicate that upvelling radiance is greater from 
locations in the dump site area which contained less then 2 milligrams per 
liter total suspended material than from non-dump site areas which contained 
up to 9 milligrams per liter. 
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INTRODUCTION 

An objective in the NASA Langley Research Center remote sensing program 
is to demonstrate that information derived from remote sensing of 
chlorophyll-ji and suspended sediment can be used to control biodegradable 
water pollution. In reference 1 it is shown that chlorophyll- a produced by 
blue-green phytoplankton can be remotely mapped with a near-infrared 
detector when the concentration of chlorophyll- is greater than 3^ tg/1 - a 
condition commonly referred to as a phytoplankton "bloom." However, as 
pointed out in reference 2, a "bloom" is a critical biodegradable condition 

r 

where additional energy input above the "bloom" level becomes a stress , and 
the system becomes poisoned. Therefore, it is desirable to detect and 
monitor phytoplankton growth below "bloom" concentrations in order to 
j/rovide biodegradable pollution control information. In reference 3 it is 
shown that below "bloom” concentrations the backscattered radiance through 
the optical filters flown was not proportional to chlorophyll- a_ concentrations , 
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and an optical filtering system for identifying the chlorophyll-ja tack- 
scattering from the total backscattered radiance is proposed. The major 
objective of the present reported photographic mission was to determine the 
effectiveness of the optical filtering system proposed in reference 3 for 
separating chlorophyll-a-containing phytoplankton radiance from total 
radiance when the concentration of phytoplankton is below f, bloom ,f conditions. 

EXPERIMENTAL METHOD 

Two photographic missions were flown at an altitude of 5.3 kilometers (km) 
by an IIA3A, Wallops Flight Center, C-5U aircraft. One on April 10, 1975, 
between the hours of 10 a.m. and 11:25 a.m. eastern daylight time and the other 
on April 13, 1975, between the hours of 9:53 a.m. and 11:03 a.m. eastern day- 
light time. Four Hasselblad cameras were used in the missions. A description 
of the camera system and the C-5J4 aircraft is included in reference k along with 
pertinent documentation of the overall remote sensing missions. The camera 
settings and film type used for the mission are shown in table 1. 

The time frame of both missions was over the later portion of the 
calculated sunglint free photographic window for a morning flight, 9 a.m. to 
11 a.m. eastern daylight time. The calculated photographic window is 
dependent upon the location of the remote sensing site, the day of the year, the 
field of view of the camera system, end the condition of water surface 
roughness. Since water surface roughness is not predictable in advance, 
calculations to determine the sunglint free window were based upon a smooth 
water surface. Any increase in sunglint effects from surface waves was 
expected to be overcome by overlapping the photographs 50 percent. In figure 
1 are shown the photographic flight lines over the New York Bight test site. 
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and the approximate location of the photographs taken during the April 10, 

1975s mission. It can be seen in figure 1 that 50 percent overlap of the 
photographs occurred on three of the five flight lines , evidently, when the 
aircraft was heading into the wind. For the other two flight lines the over- 
lap was on the order of 20 percent . ■ 

The Hasselblad cameras were equipped with different optical filters. 

The spectral transmittance (T) of each and the type used is shown in figure 2. 
Additional information concerning the filters can he found in references 5 
and 6. 

The bandwidth and location of the optical filters were selected based 
on the following: 

The narrow-band filter centered at 525 nanometers in figure 2 integrates 
the radiance around the central wavelength where, as shown in reference 7* 
the ratio of the diffused backscattering to diffused absorption coefficient 
of chlorophyll-si is equal to that for water. It is also shown (ref. 7} that 
for wavelengths less than this central wavelength the ratio of diffused 
backscattering to absorption coefficients of chlcrophyll-a is less than the 
ratio for water, and ax wavelengths greater than this central wavelength the 
ratio is greater. Thus, since the diffused backscattering to absorption 
coefficient ratio is proportional to reflectance (ref. 7) s it is expected 
that the radiance integrating narrowband filter centered at 525 nanometers 
will not detect radiance variations from chlorophyll-^ in water. Therefore, 
it is assumed that it will be useful for detecting variations in radiance 
from suspended matter, including other substances in the phytoplankton that 
contain the chlorophyll-a. In reference 7 the wavelength at which the ratio 
of diffused hackscattering to diffused absorption coefficients of a particulate 
is equal to that for water is called the hinge point wavelength. 
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In figure 3 transmission spectrums of two concentrations of blue-green 
phytoplankton Anacystis Marinus , relative to its suspension solution 9 
illustrate the concept of hinge point wavelengths. The transmission spectrums 
were obtained from a Cary lit spectrophotometer (ref. 8). The transmission 
spectrum of Anacystis Harinus for the heavier concentration is represented 
by the solid curve. It can he seen in figure 3 that the transmission spectrum 
of the heavier concentration of Anacystis Marinus phytoplankton crosses the 
transmission spectrum for the lesser concentration at three wavelengths. 

These hinge point wavelengths, illustrated in figure 3, are at approximately 
520, 570 , and TOO nanometers. At wavelengths less than 520 nanometers the 
relative transmission decreases, between the wavelengths of 520 and 570 
nanometers increases, again decreases between 570 and 700 nanometers, and, 
finally, increases again beyond 700 nanometers with increasing concentration 
of chlorophyll a, in Anacystis Marinus phytoplankton. The hinge points 
Illustrated by the transmission spectrums in figure 3 are not necessarily 
located at the same wavelengths as hinge points from reflectance spectz'ums, 
because transmission spectrums can consist of both diffused and direct light 
while reflectance spectrums consist of only diffused light. Thus, the central 
wavelength for the narrowband filter that must normalize chlorophyll 
reflectance with varying concentration was selected at 520 nanometers from 
results of the study in reference 7 and not from the transmission spectra 
shown in figure 3. This optical filter is labeled 525 B-3 filter and called 
the suspended sediment filter in figure 3. 

It is further shown in reference 7 that the backsc&ttering radiance at 
560 nanometers from chlorophyll ja concentrations between 1 and 30 lig/1 is 
always greater than the radiance from water. This is also indicated by the 


transmission spectra in figure 3, because between tbe wavelengths of 520 
and 570 nanometers transmission relative to water increases when chlorophyll a_ 
increases. The narrow-band filter with the central wavelength located at 
55^ nanometers should be able to detect the change in upwelled radiance from 
varying concentrations of chlorophyll a_ below the chlorophyll a concentration 
of 30 pg/1, However, as shown in reference 3, the upwelled radiance from the 
chlorophyll ia would he combined with the upwelled radiance from all other 
substances in the water, and must be separated from these substances in order 
to be useful for mapping chlorophyll a in a body of water. Therefore, it is 
surmised that the upwelled radiance from chlorophyll a concentrations below 
3^ Jig/1 can be obtained by a differencing or a ratio technique utilizing the 
radiance values from the two narrow-band filteres shown in. figures 2 and 3. 

The 89 B near-infrared filter in figure 2 must be flown simultaneously 
with the two narrow-band filters to determine if the water area under 
surveillance is phytoplankton '’bloomed, 11 because as shown in reference 3, the 
two narrow-band filter systems would be required for "imbloomed" phytoplankton 
water. In the "bloomed" water areas the concentrations of chlorophyll a_ could 
be determined from the near-infrared system alone (ref. l). 

The 12 yellow filter in figure 2 has two purposes. It has potential for 
identifying phytoplankton color groups in phytoplankton "bloomed" water, as 
it did in reference 3 for blue-green phytoplankton, and the 12 yellow filter, 
along with the other filter systems, has a potential for detecting radiance in 
"nonbloomed" phytoplankton water from large particulates , because it fills 
the radiance spectral gap, between 580 and TOO nanometers, not covered by the 


other filters. 
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Increasing concentrations of large particulates are expected to increase 
reflected radiance from "nonbloomed” phytoplankton -water In the yellow and 
red regions of the spectrum^ This increased reflected radiance from large 
particulates is shofm theoretic ally in reference 9 and is believed to be 
seen through the 12 yellow filter in photography over the Potomac River 
(ref. 3). 

RESULTS ABD DISCUSSION 

Presented In figures k end 5 are 276 positive prints of photographs that 
were taken with the Hasselblad cameras over the Hew York Eight from 5-3 
kilometer altitude on April 10 and 13, 1975 » The negatives of these prints 
are retained at the Chesapeake Bay Ecological Program Data Center, Wallops 
Island, Virginia, 23377- The photographs taken on April 10 are shown in 
figure 4 and those taken on April 13 in figure 5- The positive prints, 
because they are reproduced from the negatives, do not show as much detail as 
can be seen in the negatives. Positive transparencies, which show more detail 
than the positive prints, but less than the negatives, were used to obtain the 
photographic film density data used in this report. On both days only tiro of 
the cameras functioned. On both days the 89B (near-infrared) Wratteu optical 
filter system operated in addition to the 55^ (green-yellow) B-3 Baird-Atomic 
optical filter system on April 10, 1975, and the 525 (blue-green) B-3 Baird- 
Atomic optical filter system on April 13, 1975- It is believed that the very 
cold temperature (-1 6 degrees C), existing in the earners compartment of the 
aircraft at the time of the mission, is responsible for the poor performance 
of the camera system* The cameras which operated, operated Intermittently. 
This can be seen in pictures 6, 7, 8 from the 55^ (green-yellow) B-3 Baird- 
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Atomic optical filter system in figure 4 and the lack of pictures 44 through 
50 for the 525 (blue-green ) B-3 BaircL-Atomic optical filter system in 
figure 5 . Because of the poor operating performance of the camera system any 
reference to location of features seen in the pictures relative to land is 
questionable, except in photographs which include land features or known buoy 
systems. In some cases, features seen in the pictures of one flight line were 
also seen in the pictures from the adjacent flight lines due to overlap of 
the photography, see figure 1 and pictures 6„ T 9 8, 35 » 36 of figure 4. When 
this happened the location of the pictures between the two flight lines 
became more exact. 

The bright hemispherical spot originating on the border of nearly every 
picture in figures 4 and 5 is the image of the Sun reflected off the water 
surface. The Sun im age is generally largest in area in near-infrared 
photographs, because the upwelling light (radiance from the water body) is 
very weak at near— infrared wavelengths. The weak upwelling light produces 
better contrast with the surface reflected Bun image in the near— infrared 
photographs than it does at other wavelengths , because at the other wave- 
lengths the upwelling light is much greater. See pictures 2, 3, 4, and 5 of 
figure 4a and b. Whenever the image of the Sun appears in a photograph, the 
remotely sensed data obtained from that area of the photograph are not useful, 
because the direct reflection of the sun off the water surface produces 
drastic effects upon the variation in off-nadir radiance that is not detected 
at other azimuth angles. 

Many of the near-infrared photographs, figures 4a and 5a » contain white 
fuzzy lines. These lines are caused by light flashes in the camera, from the 
discharge of static electricity that builds up on 2424 near- infrared film. 


139 


A "bright area, which appears to he symmetrical about the center of e&ch 
picture, can be seen in all the scenes of figures 4b and 5b. The intensity 
of this bright area increases, in general, in photographs taken along the 
flight line as the aircraft heads toward land, and, decreases when the air- 
craft heads toward the open sea. See pictures 1 through 20, 33 through Hf, 
and 6l through 74 with the aircraft heading towards land, and 21 through 32 
and 48 through 60 with the aircraft heading out to sea in figure 4b. The 
intensity of the bright area, in general, is strongest near the center of 
each photograph and diminishes with increasing distance from the center of 
each photograph. It is shown in reference 10, from densitometer traces over 
photographs exposed from 3 kilometers altitude through two of the Hasselhlad 
cameras used in this photographic mission, that the center of photographs are 
exposed to a higher degree than the edge. This effect is normal in photo- 
graphy, and is caused by camera-lens, off-axis, falloff, which in some camera 
systems must be corrected by as much as l/(cos G) , where 0 is the angle 
off the principal axis of the lens (ref. 11), However, for wide angle remote 
sensing the camera-lens, off-axis, falloff is not .the only off-axis variable, 
because atmospheric back scattering also varies from the center of the picture 
to the edge of the border. The a+r backscattering increases from 
the center of the photograph to the edge of the film, and the increase is 
dependent upon increasing distance from the camera to the point of concern in 
the water. Since the increase in distance, similar to camera falloff, is 
dependent upon 0, the increased off-axis backscattering of sunlight compen- 
sates somewhat for the camera-lens, i.e., off-axis, falloff. Thus, the 
calibration of off-axis falloff must be performed during the photographic 
mission to account for atmospheric effects for a particular day. Whenever 
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off-axis radiance data are presented in this report it will have heen 
normalized to nadir "by multiplying it by a correction factor that is dependent 
upon the off- .adir viewing angle. The correction factors were obtained for 
each filter-film-camera system from densitometer traces through the center of 
photographs 6l {a and b) on April 10, 1975, and ik (a and b) on April 13, 

1975* For the reason previously explained, these correction factors cannot 
be applied in areas of the photographs containing the Sun’s image. 

In addition to camera-lens falloff and atmospheric back scattering, a third 
element affects the radiance values that expose the photographic film in 
combination with upwelling radiance from the water body. The third element 
is reflectance off the water surface of diffused sunlight. In reference 9 
it is shown that reflectance of diffused sunlight from a water surface is 
essentially constant with a variation in Sun altitude when the altitude cf 
the Sun is above 3^ degrees. Since the altitude of the Sun was above 3*J 
degrees for both photographic missions discussed in this report, it is 
assumed that the reflectance of diffused sunlight from the water surface is 
essentially constant in all photographs. 

Because the altitude and camera system is the same for each flight line, 
the camera-lens falloff, the atmospheric backs cattering, and the diffused 
reflectance off the water surface were assumed constant in consecutive 
photographs along the same flight line. It then follows from, these 
assumptions that the general increase in brightness at the center of each 
photograph (nadir) along the same flight line with the aircraft beading 
toward l^nd, and the decrease, when the aircraft is heading towards the open 
sea. must result from increasing or decreasing upwelling radiance. It is 
further assumed that the general decrease in upwelling radiance along each 
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flight line as the aircraft heads out to sea results from the settling of 
reflecting substances being discharged from the rivers that flow into the 
Hew fork Bight water. 

In some photographs additional light scattering features can be seen in 
the water. Shis is most evident in pictures 6, 7» 8, 35 s 36, 39» 40, 57 and 
58 of figure 4b, and pictures 9, 10 , 34 and 35 of figure 5b. These features 
and their approximate locations are more easily seen in figures 6, 7 and 8. 
Figures 6, 7 and 8, along with figure 9, are composite figures made from 
their corresponding photographs shown in figure 4 and 5* Figure 6, which is a 
composite made from the photographs in figure 4b, shows the features best. 
However, the features seen within the circle labeled sludge dump area in 
figure 6 are seen, because their reflectance is less than, while the features 
in the area labeled acid dump area are seen , because their reflectance is 
greater than the background water. In figure 7s a composite made from the 
near infrared photography of figure 4a, the same features can also be seen. 
However, the features within the circle in figure 7 are now seen, because 
their reflectance is greater than the background water, while in figure 6 
the same features were seen, because their reflectance was less than the 
background water. In figure 7 the features in the acid dump area are seen, 
because their reflectance is slightly greater than, while in figure 6 these 
same features have a much greater reflectance than the baclcground water. 
Therefore, the substance being dumped within the circle in figures 6 and 7 
must he a different substance than that in the acid dump area, because it 
has different reflecting properties 2 Native to the background water through 
the two broadband optical filters. Since the location of the circled area 
in figures 6 and 7 relative to the area labeled acid dump area is consistent 
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with the designated dump areas for waste sludge and waste acid, shown in 
reference if, it is assumed that the substance being dumped within the circle 
in figures 6 and 7 is waste sludge and the substance in the acid dump area 
must be acid residue. If this assumption is correct, it should be possible 
to distinguish waste sludge from waste acid by a remote sensing technique. 

It is also apparent that the substance being dumped within the circled area 
of figures 6 and 7 has a hinge point wavelength, because at some wavelength, 
between the broadband optical filters used in this photography, the substance 
must have the same reflectance properties as the background water. The 
significance of a hinge point wavelength for the separation of one substance 
in water relative to another using a remote sensing technique is discussed 
in reference 7. 

In figure 8 brightly reflecting features are visible from the water 
through the 525 (blue-green) B-3 Baird-Atomic optical filter. These bright 
features are labeled as drifting acid residue in lieu of sludge residue, 
because their location coincides more with the acid dump area than, it does 
with the sludge dump area, and their reflectance is, as it was in figure 6 
through the 55^ (green) B-3 Baird-Atomic optical filter, much greater than 
the background water. The fact that its reflectance is greater than the 
background water through the 525 ("blue-green) B-3 Baird-Atomic filter is 
important, because it relates to the assumed signature of acid and sludge 
obtained from figures 6 and 7» and illustrates the significance of the 
spectral response of reflecting substance for specific identification of 
these substances. 

The sludge dump area cannot be identified in the photographs taken on 
April 13 s 1975, because there are no easily identified areas in figure 8 where 
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visual observation shows reflectance from substances to he less than, and in 
figure 9 greater than, the background water. 

In figure 10 is presented the variation of nadir radiance through the 

55^- (green) and 8pB (HUt) optical filters along Xbran line 4640 on April 10, 

/ 

t 

19T5, versus miles from Staten Island. The radiance is in the form of 
relative film transmittance, obtained from densitometer traces over the 
center of positive transparencies of photographs 33 to 4?, figures 1 and 4, 
from which has been subtracted relative film transmittance through the 
unexposed portion of the film. The nadir radiance values are unique in 
remote sensing, because they result from the perpendicular viewing distances 
through the atmosphere and into the water. As such, relative to each other, 
they require no off-nadir atmospheric or cairera-lens corrections. In addition, 
since they are always in the plane of the sun, any effect upon radiance values 
from the variation in sun azimuth angle should be normalized. 

In figure 10 the value of at each location through the 554 (green) 

optical filter is always greater than the value of R through the 89B (NIB ) 
optical filter. Through each optical filter the value of is greatest 

over land, and least for the two data points beyond 33 nautical miles from 
Staten Island, labeled uniform water in figure 10. The two constant values 
of R n through each optical filter beyond station 33 were obtained from 
photographs 33 and 34 in figure 4. The dashed lines, drawn through the 
minimum, constant values of R , in figure 4, represent the radiance values 
for uniform water. Thus, the dashed lines in figure 4 are base reference 
values for increased upwelling radiance resulting from suspended materials 
during this remote sensing mission. 
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In figure 10 there is a general decay in R^ through hoth optical filters 
as distance from land increases- The general decay in R^ is distributed over 
two areas; one area where the photographs showed that heavy clouds existed, 
and the other is over the waste dump area. The increase in R^ over the 
waste dump area, between 22 and 3^ nautical miles from Staten Island, is much 
more prominent through the green filter system than through the near-infrared 
filter system, but the Increase in R^ over the cloudy area, "between 12 and 
15 nautical miles from Staten Island, is more prominent at near-infrared wave- 
lengths. Thus, there is a spectral difference in reflectance between clouds 
and the substances being dumped in the waste dump areas. 

The increase in R over the waste dump area indicates, from the stand- 
n 

point of biodegradable pollution, that a substance being dumped in the waste 
dump area is reflecting sunlight out of the water that , otherwise , might be 
available for photosynthetic activity. The extent of the effect of the 
waste dump area upon reflected sunlight can be seen better In figure 11. 

In figure 11 R^ from each over-water photograph obtained through the 
55^ (green) optical filter on April 10, 19T5, is plotted, for each loran 
line, against distance from Staten Island. Also drawn on the figure are two 
dashed circles. The smaller of the two circles encloses the area where 
visible dumping is taking place during the photographic mission, and the 
larger circle is the area where the radiance from features in the water is 
very prominent in figure 6 * The location of the smaller circle is also 
consistent with the location of the designated waste sludge dump area, and 
the larger circle with the designated acid waste dump area. 

In figure 11 there is a high value of K n along each loran line that 
increases in value from loran line ^55C to loran line ^6l0. At loran line U6lC 
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R n peaks, and then decreases again toward loran line 4670. This ridged value 

of E is shown by the solid line across the data. The peak value of E on 

n n 

loran line 46lO appears to he inside the designated acid waste dump area. 

In front of the high value of R on each loran line there is a 

n 

depressed value of R^, The depressed value of R^ decreases in value from 
loran line 4550 to loran line 464o and then increases again to loran line 4670. 
The depressed value of R^ is shown by the dashed line across the data. It 
appears from figure 11 that the minimum depressed value of R^, alcng loran 
line 4640, is in the designated waste sludge dump area. It also appears, 
since the depressed value of R^ precedes the increased value of R^ along 
each loran line, that at optical wavelengths transmitted "by the 55^- (green) 
filter the upwelled radiance from the sludge waste area is less than from the 
acid waste area. This is the same conclusion arrived at from analysis of 
figures 6 and 7, "but, in figure 11, it appears the influence of the dumping 
of sludge and acid waste upon R^ is much more extensive in area then the 
designated dump sites. It is also seen in figure 11 that the photographic 
mission did not overfly the complete area influenced by waste dumping, since 
the two outer flight lines still experience increased values of at 

distances out to sea consistent with the location of the waste dump areas. 

The effect of waste dumping upon the variation of E q did not change 
appreciably from April 10, 1975 , to April 13, 1975, because, as shown in 
figure 12, the variation in R n along each loran line obtained from the 525 
(blue-green) optical filter system shows the depressed value of 
preceding an increasing value of R^ when the waste dump areas are overflown. 
R^ from each over-water photograph through the 525 (blue-green) optical 
filter on April 13, 1975, is shown in figure 12 by the square symbols, and 


R^ from each over-water photograph through the 554 (green) optical filter 

on April 10, 1975a is represented "by the circles. R n is plotted against 

miles from Staten Island in a staggered formation to isolate R^ .ong 

each loran line from both days. The dashed lines in figure 12, shown and 

discussed for figure 10, are the minimum values of R for each filter-film 

n 

camera system. The minimum value of was different for each filter-film 

camera system, and on April 13, -1975s "was obtained on only loran flight 
lines 4614-0 end 4670, perhaps, because the flight lines did not extend as far 
out to sea on that day as they did on April 10, 1975- Nevertheless, it is 
clear in figure 12 that the influence of waste dumping upon was as 

extensive on April 13 as it was on April 10. 

Since not a single in situ data area visited by the boat or the heli- 
copter occurred in the center of a photograph, it was impossible to relate 
H n to in situ data in this experiment. This was expected prior to the 
photographic mission, because in an area as extensive as the New York Eight, 
a photographic mission from 5.3 km altitude cannot be preplanned .to pinpoint 
the center of a photograph to a specific location. In future experiments , 
the calibration of the optical filter system with in situ water data might 
better be accomplished from a helicopter. Prom a helicopter the remote 
sensing data could be taken over a location at the same time as the in situ 
water data. 

Even off -nadir radiance data, R, was not obtained in this experiment 
from many of the in situ data point locations, because the lack of 50 percent 
overlap in the photography placed the boat , and some of the helicopter in situ 
data point locations in the sunglint areas of the photographs. In all, only 
eight in situ data points were considered useful for analysis of the 
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photographic radiance data. These -were from the helicopter, and were 

restricted to April 13, because the helicopter did not fly on April 10. On 

April 13 two types of in situ data were obtained from helicopter dip samples; 

the concentration of chlorophyll- a_ and the total suspended material weight. 

% 

Since two of the photographic systems did not operate on April 13, the 
calibration of the optical filtering system could not be accomplished with 
the chlorophyll-a_ data. 

Presented in figure 13 is the variation of E through the 525 (blue- 
green) optical filter obtained from the in situ data point locations on 
April 13, 1975- R is plotted versus measured total suspended material 
weight in surface samples obtained from these same locations. Data points 
from the waste disposal areas are represented by the circles, and those 
considered to be outside of the waste disposal areas by the squares. 

It is seen in figure 13 that the variation of R in this experiment is 
not very sensitive to total suspended material weight measured in surface 
samples, because with a factor of 10 variation in total suspended material 
weight B did not vary appreciably. It is also apparent that the remote 
sensing data shown here is affected more by other parameters in the water, 
since R is greatest in the waste disposal areas where the total suspended 
material weight in surface samples was ihe least, 

CONCLUDING EEMABICS 

The major objective of this photographic experiment, the remote, 
quantification of chlorophyll- a_ in chlorophyll -a_ containing phytoplankton 
when phytoplankton "bloom" conditions do not exist, was not accomplished, 
because of camera system malfunctions. However, some important results, 
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concerning the New York Bight waste dump area, were obtained from analysis 
of radiance data through the camera systems that did operate. 

Comparison of broadband black and white photography , exposed through an 
optical filter that transmits only near-infrared wavelengths (TOO to 900 nm) 
with photography exposed through an optical filter that transmits only green- 
yellow wavelengths (530 to 580 nm), shows that a substance being dumped in 
the sludge dump area, assumed to be sanitation sludge, has a different 
spectral signature relative to the background water than does a reflecting 
substance from the acid dump area, assumed to be acid waste residue. 

The substance being dumped in the sludge dump area has a greater 
reflectance at near-infrared wavelengths and a lesser reflectance at green- 
yellow wavelengths than i e background water. 

The substance in the acid dump area has a greater reflectance than the 
background water through both optical filters. In addition, the reflectance 
from the substance in the acid dump area, is much greater in the green-yellow 
spectral region than it is in the near-infrared region. 

Film density data, obtained from the green-yellow wavelength photography , 
show that waste dumping in the Hew York Bight increases the backscattered 
sunlight from the water over a much greater erea than the designated dump site 
areas. It is clear from the photography that acid waste is the primary 
substance responsible for the increased backscattered sunlight , but the 
green-yellow wavelength data suggests that sludge dumping is also 
contributing. 

Film density data, from locations on photographs where in situ measure- 
ments of total suspended material weight in surface samples were obtained, 
indicate that upwelling radiance from areas in the dump site area containing 




less than 2 milligrams per liter total suspended material weight is greater 
than radiance from non-dump site areas that contained values of total sus- 
pended material -weight up to 9 milligrams per liter. 
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TABLE 1.- PHOTOGRAPH SENSOR COTOPLEMENT Alii' CAMERA SETTINGS 


Focal Film 


Camera 

length 

mm 

(in) 


Filter 

format 

mm 

(in) 

Film type 

Speed 

(sec) 

f 

timber 

Hasselblad 

40 

(1.57) 

554 

{ green-yellow ) a 

70 

(2.76) 

2402 Black & White 

1/250 

5.6 

Hasselblad 
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Figure 2 .- The spectral variation of transmission, T, for the 525 central wavelength, B-3 Baird-Atomic, the 
55^ central wavelength, B-3 Baird-Atomic , the 12 Wratten, and the 89 B Wratten optical filters. 
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Figure 3.- Transmission spectrums of tvo concentrations of "blue-green phytoplankton 
Anacystis Marinus , relative to its suspension solution, the transmission of the 
525 and 55 ^ B-3 Baird-Atomic optical filter, and the transmission of the 12 
and 89 B Wratten optical filter. 
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Figure 10,- Variation of nadir radiance through the 554 (green-yellow) and 89 B ( non-inf rared) optical 
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Figure 11.- Variation of nadir radiance from all photographs taken over water through 
the 55^ (green-yellow) optical filter on April 10, 1975. 
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RABSCAT Wind Measurements Over 
The I Jew 1'ork Bight 
By 

L. C. Schroeder, W. L. Jones, Jr., end J. L Mitchell 
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RADSCAT WIND MEASUREMENTS 
OVER THE NEW YORK BIGHT 

By 

L. C. Schroeder, W. L. Jones, and J. L. Mitchell 
INTRODUCTION 

On April IT* 1975 > the Advanced Applications Flight Experiment (AAFE) 
ra&iometer/scatterometer (RADSCAT) experiment was conducted in the New York 
Bight as a part of the NOAA/NASA MESA Program* The objective of this 
experiment was to demonstrate the viability of the radar remote sensing 
technique for measurement of ocean surface wind vectors and to provide these 
measurements at selected sites for the investigation of wind-induced ocean 
currents in the bight. This report describes the RADSCAT measurement 
technique, presents the inferred wind vectors, and compares the radar derived 
wind vector with "surface truth" (tower wind vector measurements). The wind- 
induced ocean currents portion of this investigation is a NOAA responsibility 
and is not discussed in this report. 

EXPERIMENTAL TECHNIQUE 

Aircraft radar scattering measurements were obtained using the AAFE 
RADSCAT 13-9 GH scatterometer* operating on the NASA Johnson Space Center's 
C-130 aircraft (NASA-929). RADSCAT was mounted to the ramp (lower cargo 
door) and was configured so that the antenna scanned across the aircraft 

♦The reader is referred to reference 1 for an in-depth discussion of the 
AAFE RADSCAT instrument. 
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ground track (fig. 1} . During the ocean measurements , the aircraft performed 
a series of counter-clockwise 3^0 degree (360°) turns (20° hank angle), while 
the antenna was pointed to the left (inward) side of the aircraft at an angle of 
50° from the aircraft vertical axis. This maneuver permitted the antenna to 
illuminate approximately the same surface areas while undergoing continuous 
azimuth rotation as illustrated in figure 2. The polarisation was transmitted 
■horizontal, received horizontal and the incidence angle was 30°* 

In making the scatterometer measurements, the quantity of interest is the 
scattering coefficient 0° which is a measure of the ocean reflectivity. For 
incidence angles greater than 25 degrees off the nadir, the mechanism for 
hacks catter is a resonant (Bragg) scattering from ocean waves whose lengths are 
of the order of the microwave wavelength. At 13.9 GIT the ocean 0° is 
directly : proportional to the amplitude of centimeter or capillary ocean waves. 
Since the generating force for capillary wave growth is the wind stress at the 
air/sea interface and since the time constant for this wave generation is a few 
seconds , these radar scatterers are in equilibrium with the local surface 
winds. Thus the radar remote sensing technique is to infer vinclspped from a 
measurement of o' 0 . 

In this experiment , radar data taken during four circle flight patterns 
at each preselected site were processed to produce the average scattering 
coefficient as a function of flight direction (azimuth), as shown in figure 3- 
Previous research (refs. 2 and 3) has shown that both windspeed and 
direction can he inferred from these anisotropic scattering signatures. As 
previously stated, the cr° is proportional to the surface frictional 
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windspeed; However, this is a difficult quantity, to measure* Historically the 

U. 

neutral stability windspeed at 19*5 meters (m) altitude has been used by 
researchers as a geophysical parameter to be correlated with radar scattering 
measurements. An example of these radar scattering signatures for three 
different ocean 19*5 m windspeeds is shown in figure 4’, where the abscissa 
is the radar azimuth relative to the crosswind direction, Note that the 0° 
monotonically increases with windspeed and each plot is a quasi -sine-wave Df 
twice the frequency of the azimuth angle. The peaks occur in the upwind and 
downwind directions and the minima in the crosswind direction; additionally, 
the upwind peak is slightly greater than the downwind peak. For the New York 
Bight RADSCAT measurements, wind vectors were inferred, using "calibration" 
curves similar to those of figure 4. During the experiment, windspeed, 
wind direction, air temperature, and sea temperature were recorded every half 
hour at the United states Coast Guard Ambrose Tower located in the New York 
Bight. These measurements were used to provide surface truth for verifying the 
RADSCAT inferred wind vectors. The mission was flown so that the first and 
last RAJDSCAT circles would occur at Ambrose Tower to provide a direct compari- 
son of RADSCAT measurements with truth data* In addition, . I’TGAA investigators 
also obtained simultaneous current measurements from buoys at four sites 
within che test area* 

MISSION RESULTS 

The mission was flown on April 17, 1975* Figure 5 shows the location of 
. sites where RADSCAT circle patterns were flown , along with a key which 

*The 19*5 meter neutral stability windspeed is the surface frictional wind- 
speed extrapolated hr* 19*5 meters using a logarithmic profile with an air-sea 
temperature differential of zero (ref. 4) . ' 
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correlates these sites with those of the mission plan. Time would not allow 
taking data at every site originally planned. 

Table I shows the wind vectors determined from analysis of the RADSCAT 

o o 

circle data. Windspeed was determined using both a mean and cr down- 
wind. The results obtained were consistent. The cr° mean values were 
obtained by averaging the windspeeds over all azimuth angles. These wind 

vector data have also been plotted in figure 6. On this figure, the sequence 
of data taking is indicated by dashed lines, and vectors indicate the magni- 
tude and direction of winds. From these, results , it can be seen that the 
wind is nearly constant until 12:12 p.m. e.d.t. From that time and site on, 
the wind diminishes in speed and the direction becomes variable. 

The wind end temperature measurements from Ambrose Tower are presented 
in figures 7 and 8, Also shown are the RADSCAT inferred wind vectors. The 
first and last RADSCAT measurements were taken at the Ambrose site. 

Accordingly, all other RADSCAT results should be considered in a qualitative 
sense. Concurrent with the first RADSCAT measurement, thr Ambrose windspeed 
(measured at UO meters altitude) dropped abruptly from l6 knots to approximately 
10 knots while the air/sea temperature difference increased from + 2°C to 
* 6°C (fig. T) Thereafter the Ambrose windspeed was constant while the 
air/sea temperature difference steadily increased to + 7°C by the completion 
of the 1 light. The large axr/sea temperature difference will result in a 
neutral stability windspeed at Ambrose Tower which is significantly reduced 
from measured windspeeds . Front reference H, the Ambrose 10-knot windspeed 
and approximately + J°C air /sea temperature difference measured at 1300 hours 
would produce a 19-5 m neutral stability windspeed of about U. 5 . knots. Thus 
both RADSCAT windspeed measurements at Ambrose Tower are in good agreement 

i 
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with the "surface truth". Furthermore, the comparison of 'wind direction shown 
in figure 8 is also in good agreement. 

CONCLUDING- REMARKS 

RADSCAT scattering data have "been used to infer windspeed and direction 
at preselected sites in the New York Bight with the intent of: 

1. Demonstrating the HADSCAT's capability as a wind sensor. 

2. Providing ground truth for current measurements made 
contemporarily with these measurements. 

Through the use of "surface truth" wind measurements at Ambrose Tower, the 
demonstration of the RADSCAT instrument was in good agreement for both 
comparisons. 

The analysis of ocean current effects is a NQAA responsibility and is 
not available at this time. 
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TABLE I. MESA WIND SPEED DATA FROM RADSCAT MEASUREMENTS (4-17-75) 


SITE 

START 

TIME 

(EOT) 

LATITUDE 

LONGITUDE 

WIND 

DIRECTION 
(OUT OF) 

WIND SPEED 
FROM 

MEAN cr 0 

(KNOTS) 

FROM 

DOWNWIND cr° 

i 

10 08 

40.45 

-73.88 

310° 

10.8 

10.5 

7 

10 35 

4a 38 

-B. 50 

330° 

10.7 
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1048 

40.20 

-a 46 
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fX 

9.9 
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1105 

40.10 

-a 56 

310° 

11.9 
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13 

II 20 

39.84 
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40.01 

-a 48 

320° 

II. 1 
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li 

II 54 

40.04 

-a 20 

305° 

8.5 


10 

12 12 

40.22 

-73. 24 

290° 

6.8 


9 

12 29 

40.40 

-a io 

245° 

5.5 

5.4 

8 

12 47 

40. 52 

-a 39 

300° 

6.3 


1 

13 06 

40.44 

-73. 90 

285° 

5.0 
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Sgure 1. RADSCAT installed on the NASA 929 aircraft 
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Hgure 4. Typical scattering data from RADSCAT circle flights. 
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RAOSCAT measurements of MESA wind vector data 
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